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TFANSLATORS'  PREFACE  TO  THE  ENGLISH  EDITION 


Within  the  past  few  years,  interest  in  the  dynamic  stability  of  elastic 
systems  has  increased.  This  interest  is  reflected  in  the  appearance  of  a 
large  number  of  papers  on  this  subject.  Unfortunately,  few  are  in  English, 
the  majority  of  these  papers  being  written  either  in  Russian  or  German. 

In  addition,  there  is  no  English  text  available  which  presents  the  mathemat* 
ical  theory  of  the  subject  with  its  applications.  However,  such  a  text  is 
available  in  Russian.  This  work,  "Dynamic  Stability  of  Elastic  Systems," 
(Gostekhizdat,  Moscow,  1956)  by  V,  V.  Bolotin,  together  with  its  German 
translation  (Veb  Deutscher  Verlag  Der  Wissenschaften,  Berlin,  1961),  is  the 
only  comprehensive  book  on  the  subject  of  dynamic  stability.  To  fill  this  need 
in  American  scientific  literature,  it  was  decided  to  undertake  the  translation 
of  this  unique  book. 

At  the  time  the  translation  was  initiated,  the  German  edition  was  not 
yet  available,  and  the  book  was  translated  from  the  Russian.  With  the 
appearance  of  the  German  edition,  the  translation  was  checked  and  all 
changes  in  the  German  edition  were  incorporated  into  the  English  translation. 
Footnotes  in  the  Russian  edition  were  eliminated  and  the  German  referencing 
system  was  adapted  for  our  use.  A  number  of  notes  were  made  by  the  trans¬ 
lators  in  order  to  bring  the  American  edition  of  the  text  up  to  date.  Refer¬ 
ences  to  new  works  published  in  the  United  States  and  an  additional  Russian 
reference  were  added. 

The  translators  would  also  like  to  thank  Dr.  Paul  Seide,  Dr.  Robert 
Cooper,  and  Dr.  John  Yao  of  Aerospace  Corporation  for  reading  the 
translation. 


El  Segundo,  California,  November  196Z  V.I.  W.,  K.  N.T.,  K.  D.G. 
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PREFACE  TO  THE  GERMAN  EDITION 


Four  years  have  passed  since  the  appearance  of  the  Russian  edition 
of  the  present  book.  During  this  time  a  series  of  works  on  dynamic  stability 
were  published,  containing  interesting  results.  In  addition,  the  general 
theory  was  applied  to  a  new  class  of  problems  whose  analysis  appeared  to 
be  of  a  mathematical  nature  which  fell  within  the  narrow  limits  of  the  theory 
of  dynamic  stability.  All  of  the  above  are  considered  in  the  German  edition 
of  the  book.  References  are  included  of  some  works  published  thru  1956, 
which  were  unknown  to  me  when  the  book  was  first  published. 

G.  Schmidt,  from  the  Institute  for  Applied  Mathematics  and  Mechanics, 
German  Academy  of  Sciences,  Berlin,  prepared  the  excellent  translation 
and  contributed  corrections  in  some  places.  The  translation  was  also 
reviewed  at  the  Institute  for  Vibration  Technology,  Karlsruhe  Technical 
University;  F.  Weidenhammer  and  G.  Bens  especially  gave  worthwhile 
advice.  C.  W.  Mishenkov,  from  the  Moscow  Power  Institute,  assisted  me 
with  the  preparation  of  complete  references.  To  all  these  people,  I  would 
like  to  express  my  deep  thanks. 


December  1960 
Moscow 


Author  -  V.  V.  Bolotin 


PREFACE  TO  THE  RUSSIAN  EDITION 


This  book  is  an  attempt  to  present  systematically  the  general  theory 
of  dynamic  stability  of  elastic  systems  and  its  numerous  applications. 
Investigations  of  the  author  are  used  as  the  basis  for  the  book,  part  of  which 
was  published  previously  in  the  form  of  separate  articles.  The  author's 
method  of  presentation  is  retained  where  the  problems  treated  have  been 
analyzed  by  other  authors. 

The  book  is  devoted  to  the  solution  of  technical  problems.  As  in  every 
other  engineering  (or  physics)  investigation,  the  presentation  consists  of 
first  choosing  an  initial  scheme  or  pattern,  and  then  using  the  approximate 
mathematical  methods  to  obtain  readily  understood  results.  This  intent, 
and  the  desire  to  make  the  book  easily  understood  by  a  large  number  of 
readers,  is  reflected  in  the  arrangement  and  structure  of  the  book. 

The  book  consists  of  three  parts.  PART  I  is  concerned  with  the 
simplest  problems  of  dynamic  stability  which  do  not  require  complicated 
mathematical  methods  for  their  solutions.  By  using  these  problems,  the 
author  wishes  to  acquaint  the  reader  with  previously  investigated  problems. 
At  the  same  time,  certain  peculiarities  of  the  phenomena  of  instability  are 
clarified,  which  previously  have  been  only  sketchily  mentioned.  PART  I 
also  contains  methods  of  solution  of  the  general  problem. 

PART  II  begins  with  two  chapters  containing  the  minimum  necessary 
mathematical  information;  a  conversant  reader  can  disregard  these  chapters. 
The  properties  of  the  general  equations  of  dynamic  stability  are  then 
examined;  methods  are  presented  for  the  determination  of  the  boundaries  of 
the  regions  of  instability  and  the  amplitudes  of  parametrically  excited 
vibrations  for  the  general  case. 
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PART  III  is  concerned  with  applications.  Various  problems  of  the 
dynamic  instability  of  straight  rods,  arches,  beams,  statically  indeterminant 
rod  systems,  plates,  and  shells  are  examined.  The  choice  of  examples  was 
dictated  by  the  desire  to  illustrate  the  general  methods  and  present  solutions 
to  practical  problems.  The  number  of  examples  was  limited  by  the  size  of 
the  book. 

I  would  like  to  take  this  opportunity  to  express  my  sincere  thanks  to 
A.  S.  Vol'mir  for  having  read  the  manuscript  and  for  having  given  valuable 
advice. 


January  1956  Author  -  V.  V.  Bolotin 

Moscow 
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ABSTRACT 


This  is  Volume  IT  of  four  proposed' volumes  of  the  translation  of 
V.  V.  Bolotin' s  book,  "The  Dynamic  St.ability  of  Elastic  Systems.  " 
Volume  II  contains  the  translation  of  Chapter  Three,  Chapter 
Four,  and  Chapter  Five.  Nonlinear  effects  are  discussed  in 
Chapter  Three.  Free  and  forced  vibrations  of  a  nonlinear  sys¬ 
tem  are  discussed  in  Chapter  Four. 

Amplitudes  of  vibrations  at  the  principal  parametric  resonance 
obtained  by  using  nonlinear  theory  are  discussed  in  Chapter  Five. 
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CHAPTER  THREE 


DETERMINATION  OF  NONLINEAR  TERMS 

•11.  PRELIMINARY  REMARKS 

From  linear  theory  one  expects  the  vibration  amplitudes  in  the  regions 
of  dynamic  instability  to  increase  unboundedly  with  time,  and  indeed  very 
rapidly,  i.e.,  they  increase  exponentially.  However,  this  conclusion 
contradicts  experimental  results,  which  show  that  vibrations  with  steady- 
state  amplitudes  exist  in  the  instability  region.  The  first  portion  of  the 
oscillogram  (Fig.  8)  shows  that  the  amplitude  of  the  vibrations  increases 
approximately  exponentially.  As  the  amplitudes  increase,  the  character  of 
the  vibrations  changes;  the  speed  of  the  growth  gradually  decreases  until 
vibrations  of  constant  (or  almost  constant)  amplitude  are  finally  established. 

The  forces  acting  upon  a  rod  can  be  considered  to  be  linear  functions  of 
displacements,  velocities,  and  accelerations  only  for  sufficiently  small 
deflections.  With  increasing  amplitudes,  the  influence  of  nonlinear  factors 
becomes  more  and  more  apparent,  i.  e.  ,  these  terms  limit  the  infinite 
increase  of  amplitudes  predicted  by  linear  theory.  Therefore,  linear  theory 
cannot  determine  either  the  magnitude  of  the  steady-state  amplitudes  or 
whether  or  not  the  vibrations  become  stationary.  Determination  of  the  above 
is  possible  only  on  the  basis  of  nonlinear  theory. 

The  question  arises  whether  the  presence  of  nonlinear  factors  introduces 
some  changes  in  the  distribution  of  regions  of  instability  given  by  linear  theory. 
This  is  not  the  case.  The  boundaries  of  the  regions  of  dynamic  instability  can 
be  given  accurately  by  the  linear  differential  equations. 

While  rigorous  investigation  of  this  question  will  be  postponed  until 
Chapter  XVI,  in  the  meantime  let  us  refer  to  its  analogy  with  the  problems  of 
static  stability,  for  example,  to  the  problem  of  the  buckling  of  a  straight  rod 
subjected  to  a  static  longitudinal  force.  If  one  linearizes  the  equations  of  the 
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problem,  it  is  well-known  that  one  obtains  the  correct  magnitude  of  the 
critical  force  (Euler  buckling  force);  however,  the  deflections  of  the  rod  in 
the  critical  and  postcritical  stages  cannot  be  determined.  1  For  the  determi¬ 
nation  of  these  deflections,  one  must  take  an  exact  (nonlinear)  expression  of 
the  curvature. 

The  idea  that  the  linear  approximation  is  sufficient  for  determining  the 
dynamic  stability  of  a  rod  was  first  clearly  formulated  by  Gol'denblat 
(Ref.  55).  Assertions  to  the  contrary  can  frequently  be  found.  Beliaev 
(Ref.  4)  assumed  that  the  consideration  of  the  nonlinear  curvature  must  in 
some  manner  be  reflected  on  the  regions  of  dynamic  instability. 

For  a  more  precise  determination  of  the  boundaries  of  the  regions  of 
instability  at  large  values  of  P^,  Beliaev  proposed  to  make  use  of  "differential 
equations  for  a  slightly  curved  rod".  Analogous  assumptions  have  been  made 
repeatedly  up  to  the  present. 

An  investigation  of  the  relationship  between  linear  and  nonlinear  theory 
is  of  basic  importance  to  the  problems  under  consideration.  We  shall  return 
to  this  question  in  the  second  part  of  the  book,  where  we  shall  relate  it  to  the 
results  of  Liapunov's  theory  of  the  stability  of  motion. 

•12.  BUCKLING  IN  THE  POSTCRITICAL  REGION 

In  the  present  chapter  we  shall  analyze  the  nonlinearities  on  which  the 
amplitudes  of  parametrically  excited  vibrations  depend.  We  will  begin  with 
the  simplest  problem  of  strength  of  materials  in  which  it  is  necessary  to 
consider  the  effect  of  a  nonlinearity,  i.  e.  ,  we  shall  begin  with  the  problem 
of  the  bending  of  a  straight  rod  by  a  longitudinal  force  that  exceeds  the  Euler 
buckling  value  (Fig.  15).  The  equation  for  the  buckling  (elastic  line)  has  the 
form 


e 


/• 


all. 


(3.  1) 
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where  —  is  the  exact  expression  for  the  curvature. 
P 


Take  the  distance 
measured  along  the  arc  of  the  deformed  rod  as  an  independent  variable, 
the  axial  deformation  along  the  rod  axis  is  disregarded, 
then  the  arc  length  coincides  with  the  x-coordinate  for 
the  undeformed  rod. 

Let  s  be  the  arc  length  and  ^  the  angle  between 
the  ox-axis  and  the  tangent  of  the  elastic  line  (Fig.  15). 

Taking  into  account  that  sin  ^  =  dv/ds  and  differentiat¬ 
ing  this  expression  along  the  arc  length  s,  we  obtain 


If 


rf*r  !<«■ 


Figure  15 


which  gives 


I  _  dr  <**•* .  • 

p  ”  dt  d**  etmt/’  ’ 


Since 


vmf  \'i  -  rtniy.  = 

we  obtain  the  following  expression  for  the  curvature 
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Equation  (3.  1)  for  the  buckling  of  a  rod  takes  the  form 


(3.2) 


The  nonlinear  differential  equation  obtained  can  be  reduced  to  a  form  more 
convenient  for  solution.  By  expanding  the  radical  in  the  binominal  expansion 


/rff'V 

)  '^  ir ' 

(v.) 

(3.3) 


Eq.  (3.  2)  can  be  rewritten  in  the  following  form: 


(3.4) 


The  first  term  in  Eq.  (3.  3)  corresponds  to  the  usual  linear  approxi¬ 
mation  in  strength  of  materials.  We  obtain  the  first  nonlinear  approximation 
by  retaining  two  terms  in  Eq.  (3.  3);  such  an  approximation  will  be  suitable 
for  small  deflections. 

/ 

If  the  longitudinal  force  does  not  greatly  exceed  the  Euler  buckling  value, 
then  the  elastic  curve  differs  only  slightly  from  the  form  of  the  first  charac¬ 
teristic  function  of  the  linear  problem, 


»(•)  =  / Kin  Y . 


(3.  5) 


We  shall  make  use  of  the  Galerkin  variational  method  for  finding  the 
unknown  deflection  f.  Variational  methods  recently  have  found  greater  and 
more  effective  application  in  different  divisions  of  the  applied  theory  of 
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elasticity  and  structur2d  mechanics;  we  shall  assume  that  the  reader  is 
familiar  with  these  methods. 

Following  the  method  of  Galerkin,  we  substitute  Eq.  (3,  5)  into  the 
left-hand  side  of  Eq.‘  (3.4)  and  require  the  resulting  equation  to  the  selected 
function  sin  ^  to  be  orthogonal: 
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0. 


With  the  known  integrals 


0 


we  obtain 


y  81*  ^  32 


for  the  determination  of  f.  Here,  as  before 


P* 


^KJ 

ji  • 


(3.6) 


One  of  the  possible  solutions  of  Eq.  (3.  6)  is  f  =  0;  this  solution  obviously 
corresponds  to  the  initial  (uncurved)  shape  of  the  rod.  The  non-zero  solutions 
can  be  found  from  the  condition 


We  recommend  Ref.  81  as  a  first  reading.  See  Refs.  82  and  83  for  readings 
in  English. 


(3.  7) 


I  + 


3  »•/*  . 

32  F  ^ 


P 

p7 


=  0. 


If  higher  powers  of  f  are  neglected  in  Eq.  (3.  7),  we  obtain  the  well-known 
.  2 

approximate  equation 


/  = 


(3.8) 


A  comparison  with  the  exact  solution  of  the  same  problem,  using 

elliptical  integrals,  shows  that  Eq.  (3.  8)  gives  good  results  for  f  <  0.  2i; 

consequently,  P  <  1.045Pj,,.  In  order  to  obtain  the  second  approximation, 
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one  should  retain  the  f  terms,  and  so  on. 

Note  that  this  equation  and  also  the  equations  for  the  higher  approxi¬ 
mations  can  be  obtained  as  a  special  case  from  the  results  of  Euler,  (Ref.  85). 
In  some  later  works  a  correction  to  the  second  approximation  was  calculated 
incorrectly  as  indicated  by  Nikolai  (Ref.  86). 

•13.  NONLINEAR  ELASTICITY 

1.  We  return  to  the  problem  of  the  dynamic  stability  of  a  hinged  straight 
rod  that  is  compressed  by  a  periodic  longitudinal  force.  We  shall  proceed 
from  the  equation 


r  +  2e /'  +  «*(!  -  2/1  ci>H9<)/+y(/. /',/")  =  0  (3.9) 

that  differs  from  Eq.  (2.  1)  by  the  presence  of  certain  nonlinear  functions  of 
displacements,  velocities,  and  accelerations,  «|>(f,  f ,  f").  The  determina¬ 
tion  of  this  function  constitutes  our  next  problem. 
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See  Ref.  84,  Volume  L  for  example. 
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Among  the  terms  entering  into  the  nonlinear  function  v)<  (f,  V,  f"),  it  is 
always  possible  to  single  out  such  terms  that  do  not  contain  derivatives  of 
displacements  with  respect  to  time.  These  nonlinear  terms  characterize  the 
nonlinear  elasticity  of  the  system.  Taking  this  into  account,  we  will  group 
all  the  static  nonlinear  factors  under  the  general  designation  of  nonlinear 
elasticity,  regardless  of  whether  they  are  of  a  geometrical  or  physical  origin. 

Nonlinear  elasticity  is  the  only  nonlinear  factor  in  the  problem  of  the 
bending  of  a  rod  by  a  force  exceeding  the  Euler  buckling  load.  The  governing 
equation  can  be  obtained  from  Eq.  (3.  9),  if  one  assumes  f  =  const,  and  elimi¬ 
nates  terms  explicitly  depending  on  time.  As  a  result,  we  obtain 


or 

(l +  i  »-(/)-«• 


(3.10) 


Comparing  Eq.  (3.  10)  with  Eq.  (3.  6),  we  can  conclude  that  the  nonlinear 
function  corresponding  to  the  complete  expression  for  the  curvature  has  the 
form 


r</)-  +*44^  +  ---) 


This  means  that  the  function  i|*(f)>  which  takes  into  account  the  influence 
of  the  nonlinear  curvature,  can  be  represented  in  the  form  of  a  series  that 
contains  odd  powers  of  the  deflection 
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(3.11) 


where 

.•»»  »»•  5  .T*  "*• 

y*  TT/*"  ’  y*  w  i*  . 

and  so  on. 

It  can  be  shown  that  consideration  of  certain  other 
nonlinearities  leads  also  to  nonlinear  functions  of  the  type 
in  Eq.  3.  11.  Let  us  examine,  for  example,  a  rod  with  a 
longitudinal  elastic  spring  (Fig.  16).  During  vibration, 
an  additional  longitudinal  force  arises  from  the  reaction  of 
the  spring: 


.l/»=-eir.  (3.12) 

Figure  16 

The  displacement  of  the  moving  end  of  the  rod  is  denoted 
by  w,  the  stiffness  of  the  spring  is  denoted  by  c,  and  the  axial  deformation 
of  the  rod  is  not  considered. 

Essentially  the  longitudinal  displacement  w  is  nonlinearly  related  to 
the  transverse  deflection  of  the  rod.  The  longitudinal  displacement  of  the 
moving  end  can  be  found  as  the  difference  between  the  initial  length  i  and 
the  projection  of  its  deformed  center  line 

w  =  /  -  r  coaf>d«  =  1  -  JVl- 

t  f 
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Expanding  the  radical  in  a  series 


( 


and  integrating  by  terms,  we  obtain 


(3.  13) 


The  first  term  of  this  series  is  frequently  encountered  in  the  bending  calcu 
lation  as  an  expression  for  the  deflection  of  the  moving  end  of  a  rod. 

If  we  substitute  into  Eq.  (3.  13)  the  expression 

*(•.  0  =  /(0»'n-Y- 


and  use  the  defined  integrals 


I  I 

“  Y>  J'con*  Y***  “  4*’ 


etc.  ,  we  find 


»»/*  .  a  . 


(3.  14) 


Taking  into  account  the  additional  longitudinal  force  in  Eq.  (3.  12),  the 
equation  of  the  vibrations  of  a  rod,  Eq.  (2.  1),  takes  the  form 
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r + 2t  /'+♦«»(  1  -  f  „  () 

(Compare  it  with  Eq.  (1.4)).  After  substituting  Eq.  (3.  14),  we  obtain 

r  -I  2e /'  +  «**  (l  -  I  -f  v*(/)  0 

where 


3  a*  m*  r. 


rU)  -  ■  /*  +  7,4  p-/.-- 1*  + 


Thus  a  consideration  of  a  longitudinal  elastic  spring  also  results  in  the 
nonlinear  function  of  the  type  in  Eq.  (3.  11).  We  will  limit  outselves  to  non¬ 
linear  terms  less  than  third  order  in  further  considerations  and  will  write 
the  nonlinear  function  in  the  form 


rU)  =•  yl* 


The  coefficient  of  nonlinear  elasticity 


3  . 

IS 


(3.15) 


Y  - 


8l» 


(3.16) 


where  the  second  term  accounts  for  the  effect  of  nonlinear  curvature.  It  is 
seen  from  this  formula  that,  with  accuracy  up  to  quantities  of  third  order 
magnitude,  the  influence  of  the  nonlinear  curvature  is  equivalent  to  the 
influence  of  a  longitudinal  spring  with  a  spring  stiffness 


c 


n*KJ 
•IP  ' 


1 


More  exactly,  it  is  the  coefficient  of  the  cubic  term  in  the  developed  series. 
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2.  Nonlinearities  of  the  above  type  often  arise  in  applications.  This 
nonlinearity  is  contained,  for  example,  in  rods  of  statically  indeterminate 
systems  which  remain  geometrically  unchanged  after  the  rods  are  unloaded 
(so-called  redundant  rods). 


In  fact,  during  the  bending  of  such  rods,  the  ends  approach  one 
another,  causing  additional  elastic  forces  to  appear  in  the  remaining  part  of 
the  system  tending  to  hinder  the  approach  of  the  ends.  Since  the  approach 
of  the  ends  described  by  Eq.  (3.  14)  is  nonlinearly  related  to  the  deflections 
of  the  rod,  the  elastic  forces  also  have  a  nonlinear  character. 


Nonlinear  elastic  forces  can  arise  also  in  systems  that  are  usually 
treated  as  being  statically  determinate.  For  example,  riveted  or  welded 
girders  are  statically  determinate  when  they  are  assumed  to  be  hinge- 
connected  and  statically  indeterminate  if  the  stiffness  of  the  hinges  is 
considered.  Still  more  complicated  relations  arise  in  spatial  constructions, 
e.  g.  ,  in  the  frames  of  metal  bridges.  In  determining  the  nonlinear  charac¬ 
teristics  of  any  rod  entering  into  the  composition  of  a  frame  structure,  the 
total  influence  of  the  spatial  construction  must  be  considered. 


The  presence  of  nonlinear  elasticity  reflects  on  the  carrying  capacity 

of  the  compressed  rods,  which  lose  their  stability  in  the  elastic  range.  As 

is  well-known,  for  single  rods  in  which  the  nonlinear  elasticity  is  dependent 

on  the  nonlinear  curvature,  even  an  insignificant  increase  of  critical  force 

results  in  dangerous  deformations.  Therefore,  the  critical  load  for  such  rods 

can,  in  practice,  be  considered  to  be 

4 

equal  to  the  limit  load.  If  a  rod  is  a 
part  of  a  statically  indeterminate  sys¬ 
tem,  its  deflection  in  the  post-critical 
stage  can  be  considerably  smaller.  To 
illustrate  this  concept,  let  us  consider 
a  simple  example  (Fig.  17). 


The  case  of  very  flexible  rods  where  large  deflections  occur  without  yielding 
and  where  creep  of  the  material  is  possible  will  not  be  considered. 
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Under  the  action  of  an  external  loading,  it  is  possible  for  element  AB 
to  become  unstable.  To  determine  the  coefficient  of  resistance  for  the  re¬ 
maining  part  of  the  construction,  we  compute  the  approach  of  end  points  A 
and  B  under  the  action  of  unit  forces.  (Fig.  17).  According  to  the  well-known 
equation  of  structural  mechanics. 


n  = 


v.?ll 

^  Kf  • 


For  simplicity,  the  areas  of  all  the  rods  will  be  assumed  to  be  equal. 
Computations  give 


=  2Un*« -I- 

hence  the  coefficient  of  resistance  is 

^  I_  _ _ KJl _ 

In  accordance  with  Eq.  (3.  16)  the  coefficient  of  nonlinear  elasticity  is 


n*oi*  I  it  4  ,\ 


The  expression  in  parentheses  can  be  written  in  the  form 


tfi 


+  I 


_ _ 


4-1 


where  \  is  the  fineness  ratio  of  the  rod.  Following  from  Eqs.  (3.  10)  and 
(3.  15),  the  deflection  of  the  middle  of  the  rod  is  determined  from  the  equation 


12 


.  ■  uji  V"  I 


where  N  is  the  longitudinal  force  in  the  rod,  is  its  critical  value,  and  k 
is  the  nondimensional  coefficient 


*(«.  A) 


Si* 


+ 1. 


Since  k(a,  X)  is  of  the  order  of  magnitude  of  the  fineness  ratio  X  (at  a  =  45° 
we  have  k  ^  0.  15X)  then  the  deflection  of  rod  AB  can  be  one-tenth  the  dis¬ 
placement  of  a  single  rod. 

Further  consideration  of  problems  dealing  with  frame  constructions  in 
the  postcritical  state  falls  outside  of  the  framework  of  the  present  book. 


3.  Hitherto,  we  have  restricted  our  treatment  to  the  case 
material  with  linear  properties.  Nonlinear  elastic  prop¬ 
erties  of  a  material  can  also  be  considered  with  the  aid  of 
Eq.  (3.11).  Since  for  the  majority  of  known  materials,  the 
stress  <r  and  the  strain  e  satisfy  the  inequality  (d  (r)/(de  )<0 
(Fig.  18),  nonlinearity  of  this  kind  in  general  will  be  "soft" 

5 

in  contrast  to  "hard"  nonlinearity  of  cases  considered  earlier. 


of  an  elastic 


The  nonlinearity  is  called  "soft"  if  the  quasi-elastic  coefficient  decreases 
with  displacement  and  "hard"  if  the  quasi-elastic  coefficient  increases  with 
displacement. 


For  example,  let  the  properties  of  the  material  be  described  by  the 
relation 


a  E{e  —  ^  e*) 

where  p  is  a  constant.  Since  plane  sections  remain  plane,  we  can  assume 


*  = 


9 


dt* 


and  we  find 


for  the  bending  moment  in  the  beam  cross-section,  where  B  is  a  new  constant. 
After  substitution  in  the  ordinary  differential  equation,  the  last  term  on  the 
right-hand  side  yields  a  nonlinear  function  for  \  <  0  of  the  form  in  Eq.  (3.  15). 

The  question  arises  whether  or  not  the  nonlinear  function  in  Eq.  (3.  11) 
can  include  all  the  possible  cases  of  nonlinear  elasticity.  Even  in  the  case 
where  the  magnitude  of  elastic  forces  depends  on  the  sign  of  the  deflections, 
we  must  introduce  in  Eq.  (3.  11)  terms  containing  even  powers  of  the  deflec¬ 
tion  (nonsymmetrical  power  characteristic). 

I  Furthermore,  nonlinear  elastic  forces  exist 

'  .  whose  dependence  on  the  deflections  cannot 

_ ^  "^^e  represented  in  the  form  of  converging 

I  power  series.  The  simplest  example  of  this 

i  I 

!  sort  is  the  rod  with  deflection  restrictors 

which  are  springs  (Fig.  19). 

Figure  19  Assume  that  the  stiffness  of  the 

restrictors  is  sufficiently  small  so  that  one 
can  neglect  the  influence  of  the  additional  support  on  the  form  of  an  elastic 


t 
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curve  (only  in  this  case  can  the  nonlinearity  be  condisered  small).  The 
nonlinear  function  then  takes  the  form 


y(/)  = 


j  0  for 


where  c  is  the  spring  constant  of  the  restrictor.  An  expression  of  this  kind 
obviously  cannot  be  represented  in  the  form  of  a  power  series. 

However,  problems  of  this  type  are  seldom  encountered,  while  at  the 
same  time,  functions  such  as  Eq.  (3.  11)  encompass  the  majority  of  problems. 
In  all  following  considerations,  we  will  confine  ourselves  to  Eq.  (3. 11). 

•14.  NONUNEAR  INERTIA^ 


Previously  we  have  considered  nonlinear  factors  of  static  origin.  In 
dynamic  problems  one  must  also  consider  nonlinear  inertia  forces  and  non¬ 
linear  damping. 


As  an  example  of  nonlinear  inertia  we  will  investigate  a  simple  problem 
and  assume  that  on  the  moving  end  of  the  rod  there  is  a  concentrated  mass 
(^ig*  20).  In  this  case  the  additional  longitudinal  force 


AP  -  w”. 


arises  during  vibration  of  the  rod.  As  before,  w(t)  denotes  the 
longitudinal  displacement  of  the  moving  end 


V 


41  '  ^  M 


I. 

b-r 

1 

a 

J 

Figure  ZO 

This  designates  the  nonlinear  (more  exactly,  cubic)  part  of 
the  expression  for  the  inertia  effect  that  enters  mechanical  problems  through 
the  choice  of  variables.  In  Lagrangian  mechanics,  on  the  other  hand,  one 
speaks  only  generally  of  inertia  forces  even  if  these  are  described  by  a  non¬ 
linear  expression  relating  the  chosen  generalized  coordinates. 
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By  considering  this  longitudinal  force,  the  differential  equation  of  the 
tranverse  vibrations  takes  the  form 


r -t- 2«/' +  a»*(l - 


p* 


0 


where 


•*"  =  [//" + im  +  T„-  -J- /*[//"  +  »(/')*]  +  •  •  •. 

If  we  substitute  the  frequency  Q  of  the  natural  vibrations  of  a  loaded  rod  and 
the  exciting  parameter  p,  we  obtain 

/"  +  2e/'  +  -  tVr-HflO/  4-  2h/|//"  +  (/')»;=<>. 

Terms  which  are  higher  than  third  order  have  been  neglected.  The  influence 
of  the  nonlinear  inertia  forces  is  considered  through  the  function 

y(/./'.n  =  2»»/r/r  +  (/Tj  (3.  i?) 


The  coefficient 


(3.18) 

will  be  called,  the  coefficient  of  nonlinear  inertia. 

2.  A  nonlinear  expression  of  the  type  in  Eq.  (3.  17)  was  first  found  in  an 
article  by  Krylov  and  Bogoliubov  (Ref.  5)  in  connection  with  the  investigation 
of  the  natural  vibrations  of  struts.  The  necessity  to  account  for  this  non¬ 
linearity  in  the  case  of  parametric  resonce  was  indicated  in  an  article  by 
Gol'denblat  (Ref.  55). 
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The  presence  of  a  concentrated  mass  on  the  moving  end  is  not  the  only 
condition  which  gives  rise  to  nonlinear  inertia  forces.  During  transverse 
vibrations,  every  section  of  the  rod  goes  through  some  longitudinal  displace¬ 
ment.  This  displacement  is  of  second  order  in  comparison  to  the  transverse 
deflection 

0 

Therefore  a  distributed  loading  due  to  the  inertia  forces 

«(*.  t)  =  — 


acts  on  the  rod.  Assume  that  this  longitudinal  loading  influences  the  shape 
of  the  vibration  only  slightly  and  again  take 


»(»./)  = 


(3.  19) 


After  substitution  and  integration,  we  find 


Obviously  the  additional  longitudinal  force  in  each  section  is 


or,  after  integration  is 

dN{x.t) 


dN{x.  t) 


I  n{f,t)d( 
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We  will  now  formulate  the  differential  equation  of  the  vibrations  of  a 
rod.  In  contrast  to  Eq.  (1. 2),  it  must  contain  an  additional  term  that 
considers  the  influence  of  the  longitudinal  force  AN(x,  t)  dependent  on  x. 
By  neglecting  damping,  the  equation  will  have  the  form 


With  the  aid  of  Gale r kin's  variational  method,  it  can  be  reduced  to  an 
ordinary  differential  equation.  In  fact,  if  we  substitute  Eq.  (3.19),  multiply 
by  sin^,  and  integrate,  we  obtain 

/»  4.  /  -f  2k/(//''  )-  (/')•!  =  0 


where 


t 


0 


The  calculation  of  the  integral  gives 


The  coefficient  x>  determined  according  to  Eq.  (3.  20),  takes  into  account  the 
influence  of  the  inertia  forces  of  the  rod  itself.  Comparing  Eqs.  (3.  18)  and 
(3.  20),  we  see  that  the  influence  of  the  inertia  forces  of  the  rod  has  the  same 
influence  as  an  equivalent  mass 


concentrated  on  the  moving  end. 
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3,  We  shall  show  how  to  estimate  the  magnitude  of  the  coefficient  of 
nonlinear  inertia  for  elements  of  complicated  rod  systems.  Consider,  for 
example,  element  AB  of  the  upper  boom  of  a  multiple  span  truss  (Fig.  21). 
Under  the  action  of  the  external  loading 

»  =  -h  ffi  e""®* 

periodic  stresses  arise  in  the  rods  of  the  truss.  The  magnitude  of  these 
stresses  can  be  calculated  by  ordinary  structural  dynamics  methods.  The 
truss  experiences  the  usual  forced  vibrations  in  which  additional  terms  of 
second  order  amplitude  arise  if  one  of  the  rods 
becomes  unstable. 

If  rod  AB  becomes  dynamically  unstable, 
its  vibrations  due  to  the  approach  of  the  ends 
(expressed  by  a  quantity  w)  will  be  accompanied 
by  vibrations  of  the  entire  truss.  Assuming  that 
the  truss  is  statically  determinate  (ideal  hinges), 
we  obtain  the  vibration  form  shown  in  Fig.  21b. 

The  amplitudes  obviously  are  quantities  of  second 
order  in  comparison  to  the  deflection  of  rod  AB. 

It  is  important  to  emphasize  that  we  are  concerned 
only  with  the  additional  deformation  of  the  truss 
arising  from  the  approach  of  the  ends  of  element 
AB. 

We  will  now  calculate  the  coefficients  of  nonlinear  elasticity.  The 
deflection  of  the  middle  of  a  span  can  be  found  from  simple  geometric 
considerations  and  is 


a) 


Figure  21 


X  = 


wL 
"4A  ’ 
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where  L  is  the  span  of  the  truss  and  h  is  its  height.  Additional  inertia 
loading  on  the  truss  as  seen  from  Fig.  21c  is 


d«(*,  0  — 


m,x 

nr  nr* 

»)  il*t» 
2*  dl*~ 


for  *5;^. 
for  *>y; 


where  denotes  the  sum  of  the  mass  of  the  truss  per  unit  length  and  the 
mass  of  the  previously  given  distributed  load.  Hence,  one  easily  finds  the 
additional  stress^  in  rod  AB. 


JN=- 


m,  Ij*  d*tt 

iHh*'  m*  • 


(3.21) 


As  seen  from  Eq.  (3.  21),  the  influence  of  the  inertia  forces  of  the 
truss  on  the  magnitude  of  the  longitudinal  force  in  the  rod  is  equivalent  to 
the  influence  of  the  concentrated  mass 

Mf 


placed  on  the  moving  end  of  the  rod.  The  coefficient  of  nonlinear  inertia  of 
the  rod  AB  is 


m,  »* 


(3.  22) 


where  m  is  the  mass  per  unit  length  of  the  rod  and  n  is  the  number  of  panels 
of  the  bottom  ginder.  The  coefficient  of  the  nonlinear  inertia  for  other  rods 


'Our  calculation  is  only  approximate;  we  actually  referred  all  the  mass  of 
the  truss  and  the  loading  to  the  lower  boom.  This  type  of  procedure  is  often 
used  in  "linear"  structural  dynamics.  It  appears  that  this  approximate 
method  gives  sufficient  accuracy  for  practical  purposes  (see  Ref.  87). 


is  determined  in  a  similar  manner.  Thus,  for  the  supporting  cross  strut 
CD  (Fig.  21),  we  obtain 


X 


Jl*(ll  —  I)*  w, 
An** m  ’ 


(3.  23) 


where  s  is  the  length  of  the  strut  and  0  is  the  angle  it  makes  with  the 
horizontal.  In  general,  if  one  considers  only  the  vertical  inertia  forces, 
the  coefficient  x  can  be  computed  by  the  following  approximate  formula  for 
beam  elements  (statically  determinate  and  also  statically  indeterminate)  of 
the  truss: 


X  — 


4mi»  f  ’ 

i 


(3.  24) 


In  Eq.  (3.  24),  4‘(s)  ordinates  of  the  influence  line  of  a  force 

which  results  from  a  unit  vertical  force  in  the  corresponding  element.  To 
derive  this  formula,  it  must  be  assumed  that  all  of  the  mass  of  the  rod 
system  can  be  referred  to  a  definite  line,  e.  g,  ,  to  the  line  of  action.  Then 
the  additional  inertia  forces  can  be  determined  by  the  formula 


J '/(.■».  t) 


,  ,  ft*W 

-W,(.v)  »)(.,) 


from  which,  according  to  the  loading  of  the  influence  line,  we  find 


f  »«,(*) 


for  the  force  in  the  rod.  The  equivalent  mass  obviously  is 


Mu  =  j  »'*(*)  ^ti¬ 


ll 


from  which  Eq.  (3.  24)  follows.  Equations  (3.  22)  and  (3.  23)  are  its  special 
cases. 

Equation  (3.  24)  can  also  be  easily  extended  to  any  rod  system;  for  this 
purpose  the  concept  of  the  "influence  line"  must  be  interpreted  in  a  broader 
sense.  Other  examples  for  the  determination  of  the  coefficient  of  nonlinear 
inertia  will  occur. 

si 5.  NONLINEAR  DAMPING 

1.  The  damping  of  free  and  forced  vibrations  has  not  yet  been  fully  con¬ 
sidered.  The  complex  and  diverse  processes  accompanying  energy 
dissipation  during  vibrations  and  the  influence  of  a  large  number  of  factors 
that  are  difficult  to  take  into  account  by  theoretical  means  have  resulted, 
up  to  this  point,  in  the  reduction  of  the  effect  of  damping  to  the  addition  of 
certain  "suitable"  terms  in  the  equation  of  the  conservative  problem.  These 
terms  are  selected  so  that  the  theoretical  results  are  in  satisfactory  agree¬ 
ment  with  experimenUil  data.  A  large  number  of  "hypotheses"  of  every  kind 
that  do  not  agree  with  the  experimental  data  are  characteristic  of  the  present 
state  of  the  art  on  damping  (Refs.  88  and  89). 

We  will  investigate  the  equation  of  damped  vibrations 

/''+«•/=«{/./')•  (3.25) 

where  R(f,  f)  is  an  additional  term  which  considers  the  resistance  forces. 
The  simplest  and  most  widely  used  expression  for  the  calculation  of  damping 
is 


«(/')= -2c/'  (3.26) 

where  c  is  a  constant,  leads  to  the  result  briefly  summarized  in  the 
following. 
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The  amplitudes  of  free  vibrations  decrease  in  a  geometric  series  so 
that  the  logarithmic  decrement  of  damping 


In -7 


/(« 


/(» +  J) 


=  «  r  =  A 


(3.27) 


is  a  quantity  that  is  independent  of  the  amplitude.  The  period  of  the  free 
vibrations  will  be  denoted  by  T: 


The  dissipation  of  energy 


AW^-jimd!  (3.28) 

for  one  period  is  proportional  to  the  square  of  the  amplitude,  and  the  relative 
dissipation 


y  = 


(3.29) 


g 

is  independent  of  the  amplitude. 

Numerous  experiments  show,  however,  that  energy  loss  during 
vibrations  is  only  slightly  dependent  on  the  velocity.  (For  example,  the  area 
of  loop  of  elastic  hysteresis  is  practically  independent  of  the  duration  of  the 
vibration  cycle,  i.e.,  of  the  speed  of  deformation.)  Meanwhile,  according 
to  Eq.  (3.26),  it  would  seem  that  the  resistance  forces  are  greatly  dependent 
on  the  velocity.  For  many  years  "the  hypothesis  of  viscous  resistance"  has 


Computations  give  =  26. 
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often  been  criticized  for  this.  It  is  proposed  here  to  consider  the  forces  of 
resistance  proportional  to  the  amplitude  and  oriented  in  the  direction 
opposite  to  that  of  the  velocity 


«(/./')-  -y|/l  •««"/': 


(3.  30) 


where  one  can  write 


•ign/'  = 


I  ' 

l-l 


for  r>0. 
for  /'  <  0, 


or  make  use  of  "the  complex  modulus."  borrowed  from  acoustics  and 
electrical  engineering,  etc. 

The  quantities  that  the  resistance  forces  depend  upon  are  not  available 
from  Eq.  (3.  26).  It  is  simply  an  indirect  method  of  considering  damping. 
Furthermore,  there  is  no  basis  at  hand  for  assigning  the  interpretation  of  a 
material  constant  to  the  damping  coefficient.  On  the  contrary,  it  is  more 
logical  (and  this  is  confirmed  by  experimental  data)  to  consider  the  corre¬ 
sponding  energy  dissipation  as  a  material  constant  or  (almost  the  same  thing) 
the  logarithmic  decrement  of  damping.  In  this  case,  the  "hypothesis  of  vis¬ 
cous  resistance"  and  also  the  methods  that  have  been  proposed  give  analogous, 
and  in  the  case  of  a  suitable  selection  of  the  constants,  coinciding  results. 

We  will  call  the  damping  linear  if  the  relative  dissipation  of  energy  tf 
is  independent  of  the  amplitude.  One  can  easily  see  that  this  definition  is 
somewhat  broader  than  the  one  that  is  generally  used.  For  example,  it  also 
includes  the  damping  described  by  Eq.  (3.  30)  which  is  generally  nonlinear. 
For  the  explicit  calculation  we  are  especially  interested  not  in  the  change  of  the 
damping  during  the  period  but  in  its  overall  effect  during  a  period.  Its  measure¬ 
ment  describes  the  relative  dissipation  of  energy.  From  this  viewpoint,  the  ex¬ 
pressions  for  the  resistance  forces  given  by  Eqs.  (3.26)  and  (3.30)  are  equivalent. 


1 


For  a  survey  of  the  literature,  see  Ref.  89. 
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The  most  detailed  investigations  on  nonlinear  damping  consider  effects 
of  the  internal  damping  in  the  material.  As  is  known,  the  energy  losses  due 
to  internal  friction  form  the  most  significant  part  of  the  general  energy  dissi¬ 
pation.  At  the  same  time,  this  part  lends  itself  to  theoretical  considerations. 

Among  the  expressions  for  the  consideration  of  internal  damping,  we 
will  indicate  only  a  generalization  of  the  "classical"  damping  expression 
Eq.  (3.  26)  by 


For  a  determination  of  constants  h  and  k,  numerous  experiments  were 
carried  out  on  bending  and  torsional  vibrations.  Thus,  according  to  the 
data  of  Lunts  (Ref.  91),  k  =  2.  17  for  steel;  according  to  other  data,  the 
index  k  takes  on  values  from  2  to  3. 

Davidenkov  (Ref.  88)  related  the  magnitude  of  internal  friction  to  the 
phenomenon  of  elastic  hysteresis.  In  order  to  describe  the  upper  and  lower 
branches  of  the  hysteresis  curve  <r  =  (r(e),  he  proposed  the  equation 

«  =  £  je  T  I  t(e.  ±  e)»  -  2- '  rj)} 

where  q  and  n  are  material  constants  and  e^  is  the 
deformation  amplitude  (Fig.  22).  On  the  basis  of  Eq. 

(3.  31),  a  series  of  special  problems  were  investigated 
by  Pisarenko  (Ref.  92). 

2.  As  is  seen  from  the  above,  the  usual  methods  for 
considering  damping  are  rarely  used  for  the  solution  of 
the  indicated  problem.  In  particular,  resulting  non¬ 
linear  "hypotheses"  are  inadequate  for  two  reasons.  First,  even  if  one 

The  physical  aspect  of  the  problem  of  internal  friction  is  discussed, 
for  example,  in  Ref.  90. 


(3.31) 
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assumes  that  they  are  correct,  they  take  into  account  only  a  comparatively 
small  part  of  the  total  dissipation  of  energy.  We  have  indicated  this  already 
(#9).  Second,  the  application  of  these  "hypotheses"  creates  serious  mathe¬ 
matical  difficulties  which  have  not  yet  been  applied  to  even  the  simplest  prob¬ 
lems  of  free  and  forced  oscillations.  All  this  stimulates  the  search  for  new 
methods  which  are  more  suita^ble  and  adaptable  to  the  calculation  of  damping. 

It  appears  to  be  more  logical  to  consider  the  coefficient  of  damping  in 
Eq.  (3.  26)  not  as  a  constant,  but  as  a  function  of  displacement.  Considering 
that  this  function  must  be  even,  and  approximating  it  by  a  power  series,  we 
obtain  the  following  expression 

-2(«  +  e,/» +  «,/♦  +  •  •  •)/'  (3.32) 

(hei'e  c,  c^,  c^i  .  •  •  constants  to  be  determined  from  experiments). 

The  expression  of  the  type  in  Eq.  (3.  32)  is  encountered  in  radio 
engineering  in  the  theory  of  vacuum  tube  generators.  This  corresponds  to 
the  so  called  "soft"  condition  of  the  generators  when  the  series  of  this 
expression  contains  the  first  two  terms,  and  to  the  "hard"  condition  (see 
Ref.  93)  when  it  contains  three  terms.  It  is  well-known  that  problems  with 
nonlinearities  of  the  type  in  Eq.  (3.  32)  are  readily  adaptable  to  mathematical 
treatment.  An  investigation  of  free  vibrations  produces  an  amplitude 
dependent  decrement  in  the  final  calculation. 

In  this  connection,  we  shall  consider  a  recent  proposal  of  Panovko 
(Ref.  89).  We  will  not  go  into  the  basic  theory  of  this  proposal  (the  elliptic 
form  of  the  loop  of  the  hysteresis  and  the  independence  of  its  area  on  the 
frequency),  but  note  only  that  it  is  equivalent  to  the  expression  for  the 
resistance  force 
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where  A  corresponds  to  the  amplitude  of  the  vibrations,  and  b  and  n  are 
constants.  We  shall  now  take 


/  =  eoafiHl  4-  i); 

where  A  is  a  constant  for  the  forced  vibration  problem  and  is  a  "slowly 
changing"  function  of  time  in  the  free  vibration  problem  with  a  small  amount 
of  damping.  Then 


Mit.n - 1) 


or 


*(/.n= (3.33) 

For  damped  vibrations  the  last  equation  is  approximate. 

In  a  certain  sense,  Eqs.  (3.  32)  and  (3.  33)  are  equivalent.  If  in 
Eq.  (3.  32)  the  damping  coefficient  is  a  function  of  the  instantaneous,  .value  of 
displacement,  then  it  is  also  dependent  on  the  amplitude.  Equation  (3.  33)  is 
simpler,  but  it  is  too  "quasi-linear".  It  becomes  unsuitable,  for  example, 
for  vibrations  with  two  different  frequencies.  Independently  from  this,  we 
are  inclined  to  favor  Eq.  (3.  32)  since  it  can  be  reduced  to  the  following 
mechanical  model. 

Let  us  return  to  our  model  (Fig.  16).  Together  with  the  usual  resist¬ 
ance  force  we  will  also  consider  a  friction  force  that  arises  during  the 
guidance  of  the  moving  support.  We  will  assume  that  it  is  proportional  to 
the  velocity  of  the  displacement  of  the  support 


* 

t 


The  equation  for  the  free  vibrations  of  a  damped  rod  takes  the  form 

r  +  +  +  (3.34) 

(Compare  with  #13  and  #14.)  By  a  termwise  differentiation  of  Eq.  (3.  14), 
we  find 

I 

f 

\ 

Substitution  into  Eq.  (3.  34)  yields 

i 

} 

+  +  +  +  + - ^0; 

which  reduces  to  a  nonlinear  function  of  the  type  in  Eq,.  (3.  32). 

If  the  displacements  are  small,  then  the  nonlinear  function  can  be 
confined  to  terms  less  than  third  order.  This  function  assumes  the  form 

y(/./') -«**/•/'  (3.35) 

whe re  c  ^  is  the  coefficient  of  nonlinear  damping 

(3.36) 


We  shall  show  that  if  a  rod  is  part  of  a  rod  system,  the  nonlinearity  of 
the  damping  can  be  explained  by  the  energy  dissipation  in  the  remaining  part 
of  the  system.  This  follows  from  considerations  similar  to  those  expressed 
in  sl4. 

As  long  as  the  amplitudes  of  the  vibrations  of  the  rod  are  sufficiently 
small,  in  practice  the  energy  dissipation  occurs  in  the  rod  itself.  Because 
of  the  approach  of  the  ends  of  the  rod,  for  larger  amplitudes  additional 
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displacements  of  the  entire  truss  that  are  nonlinearly  related  to  the 
displacements  of  the  rod  arise.  The  energy  loss  due  to  these  displacements 
also  produces  damping  which  is  nonlinear.  Under  certain  additional  assump¬ 
tions,  the  nonlinear  damping  for  this  case  can  be  determined  by  analytical 
means. 

Let  us  investigate  the  truss  shown  in  Fig.  21.  As  is  often  done,  we 
will  represent  the  resistance  forces  directed  opposite  to  the  motion  in  the 
form  of  a  distributed  loading  that  is  proportional  at  every  point  to  the  velocity 
of  the  motion.  Considering  that 


we  obtain: 


k,{L-x) 

■  2* 


w 


for 

for 


«  Si 
»  > 


1 
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where  is  the  resistance  coefficient  for  the  truss.  An  additional  longi¬ 
tudinal  force 


48*' 


w 


is  formed  in  the  element  AB.  Thus,  the  damping  in  a  truss  is  equivalent 
to  the  damping  caused  by  the  resistance  developed  in  the  moving  support 
that  has  a  coefficient 


ki  — 


48*' 


The  coefficient  of  nonlinear  damping  will  be 


‘  “  I»2m*'{> 


29 


or  approximately 


«x  •• 


J*"' 

2»**' 


(3.  37) 


where  n  is  the  number  of  panels  of  the  lower  girder.  It  is'  seen  from 
Eq.  (3.  37)  that  the  nonlinear  part  of  the  damping  increases  rapidly  as  the 
number  of  panels  of  the  truss  increases. 

3.  The  nonlinear  resistances  described  so  far  are  characterized  by  the 
fact  that,  depending  upon  the  amplitude,  its  work  in  a  period  increases  faster 
than  the  work  of  viscous  (velocity  proportional)  friction.  Somewhat  apart 
from  this,  we  have  so-called  dry  friction,  which  usually  is  considered  con¬ 
stant  in  magnitude  with  its  direction  opposite  to  that  of  the  velocity: 


«(/')= 


It  is  readily  seen  that  the  work  of  dry  friction  during  a  period 

dll' 

is  proportional  to  the  first  power  of  the  amplitude,  i.  e.  ,  it  increases  slower 
than  the  work  of  linear  viscous  friction.  Therefore,  dry  friction,  taken 
separately,  cannot  explain  the  presence  of  steady- state  vibrations  at 
parametric  resonance. 

Dry  friction  can  arise  in  support  attachments  where  it  acts  with  viscous 
resistance  (combined  friction).  It  is  interesting  that  the  presence  of  dry 
friction  in  the  moving  support  produces  damping  of  the  transverse  vibrations 
of  the  rod  according  to  a  law  that  is  characteristic  for  a  linear  resistance.  In 
other  words,  dry  friction  in  a  moving  support  and  viscous  resistance  of 
transverse  vibration  are  in  a  certain  sense  equivalent.  In  fact,  the  frictional 
force  in  the  moving  support  is 
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.\P  =  —  t,  Hiunw' 


1 
or 


l/‘=-  (//'). 

A  substitution  in  the  equation  of  free  vibrations  yields 


But  since 


=*  |/|  Hign/' 


this  equation  takes  the  fornn 


/"  +  «**/  +  -^|/|»ign/'  =  0. 


Thus,  the  presence  of  dry  friction  in  the  moving  support  is  considered 
in  the  transverse  vibration  equation  by  the  component  of  the  type  in  Eq.  (3.  30), 
i.  e.  ,  by  a  resistance  equivalent  to  the  viscous  resistance. 


TT 


Terms  of  higher  order  obviously  do  not  influence  the  sign  of  w' . 
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CHAPTER  FOUR 

J 

I  FREE  AND  FORCED  VIBRATIONS  OF  A  NONLINEAR  SYSTEM 

i 

j 

r 

j  •16.  METHOD  OF  SLOWLY  CHANGING  AMPLITUDES 

1.  The  modern  methods  for  the  investigation  of  nonlinear  vibrations  are 
based  on  the  well-known  works  of  Liapunov  (Ref.  94)  and  Pioncare'  (Ref.  95; 

;  see  also  Ref.  96). 

A  rigorous  examination  of  nonlinear  differential  equations  in  the 
general  case  leads  to  serious  mathematical  difficulties.  However,  there 
exists  a  broad  class  of  differential  equations  which  can  be  solved  by  effective 
approximate  methods.  To  this  class  of  equations  belong, in  particular,  the 
equations  describing  the  vibrations  of  systems  with  small  nonlinearities. 

Let  us  consider  the  nonlinear  differential  equation 

( 

/"  +  2e /'  +  o>» /  +  y (A /',  /")  =  0. 

By  introducing  a  nondimensional  time  and  a  nondimensional  displacement, 
we  can  write  this  equation  in  the  form 


The  prime  denotes  differentiations  with  respect  to  the  nondimensional  time. 
We  will  call  the  nonlinearity  small  if  the  condition 


n 


(4.2) 
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is  fulfilled.  In  addition,  we  will  assume  that  the  damping  is  small: 


m  n 


(4.3) 


As  we  will  see,  the  problem  under  consideration  will  satisfy  these 
conditions. 

2.  One  of  the  simplest  methods  for  the  solution  of  differential  equations 
with  small  nonlinearities  is  the  method  of  slowly  changing  amplitudes 
(van-der-Pol  method).  Although  it  is  not  sufficiently  rigorous,  it  possesses 
physical  clearness.  This  method  is  widely  applied  in  radio  engineering 
where  it  has  helped  to  obtain  a  whole  series  of  important  results  (Ref.  93). 

The  basic  idea  of  the  method  of  slowly  changing  amplitudes  is  as 
follows.  If  the  nonlinearities  and  the  damping  of  the  system  is  sufficiently 
small,  the  solution  of  the  nonlinear  equation,  Eq.  (4.  1),  at  least  during  a 
period,  differs  only  slightly  from  the  solution  of  the  linear  differential 
equation 


r  +  a)*/  =  0. 


i.  e.  ,  from  the  harmonic  vibration 

/(I)  =2  o  Bind)  i. 


It  is  said,  for  such  a  case,  that  the  solution  of  the  differential  equation, 
Eq.  (4.  1),  is  of  an  almost  periodic  character.  Accordingly,  we  will  seek  an 
approximate  solution  of  Eq.  (4.  1)  in  the  form 


/(()  =  a(<)  Binw  (; 


(4.4) 
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where  u  is  a  frequency  which  in  general  differs  from  the  "linear  frequency" 
b>,  and  a(t)  is  the  "slowly  changing  amplitude.  "  The  last  expression  can  be 
used  only  when  the  increase  in  amplitude  for  the  period  is  sufficiently  small 
in  comparison  with  its  mean  value,  i.  e.  , 


1 

^1 

w  ' 

: 

a" 

1  — cl. 

!  • 

Since  the  initial  instant  of  time  is  arbitrary,  we  will  consider  the  initial 
phase  angle  equal  to  zero. 

We  now  substitute  Eq.  (4.  4)  into  the  nonlinear  function  i|/(f,  f,  f")  and 
expand  it  in  a  Fourier  series 

y(/. /'./")  ==  <P(o,  <5)  sinw I  -f-  f(a,m)oosail  +  •  •  ’.  (4.  6) 

Terms  containing  the  higher  harmonics  are  not  wVitten  out.  In  determining 
the  coefficients 


tn 

0(a,  w)  =  J  y(/,  /")ainatdl, 
0 

in 

m 

•/'(a,  5)  =  ^  fyif,  r,  ntmaldi 


we  make  use  of  the  inequality,  Eq.  (4.  5),  which  is  approximately  written  as 

/'  ^  a  a  cnsu  t  +  a’  Bind)  (  m  a  di  coaS  I, 

/"  =  — a  d»*  sinw  t  -f  2a' d>  cosdi  t  +  o"  sindi  I  —a  at*  siiidi  I 


(4.7) 
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A  substitution  of  Eqs.  (4.  4)  and  (4.  6)  into  Eq.  (4.  1)  yields 


—  «  w*8iiia5  ( -f-  2a'(u  COSO)  (  -f-  o"  Hinw  t  2c  o'  sinw  t  +  2c  oto  co8(o< 

rt  «iiii5  {  -|-  0(n,  o>)  sind)  t  -f-  !P(o,  w)  cosS  /  +  •  •  •  ==  0. 


On  the  basis  of  Eqs.  (4.  3)  and  (4.  5),  the  underlined  terms  can  be 
neglected.  By  equating  the  coefficients  of  sin  u>t  and  cos  ut  to  zero,  we  ob¬ 
tain  the  equations 


(w*  --  <0®)o  +  0(O,  0))'=  0, 


7i 


(4.8) 


The  first  equation  of  Eq.  (4.  8)  establishes  the  dependence  of  frequency 
on  the  vibration  amplitude;  the  second  equation  determines  the  amplitude 
change  with  respect  to  time.  In  particular,  for  the  linear  system  with  damp¬ 
ing,  we  obtain:  ^ 


ci  ^  o>  ■  ■  const,  a  ---  o,  e"'*. 

3.  Another  basis  for  this  method  given  by  K.  F.  Teodorchik  (Ref.  93) 
will  now  be  discussed.  The  differential  equation,  Eq.  (4.  1),  will  be  written 
in  the  form 


r  +  «*/-(«*-  «*)  /  -  2«/'  -  y (/,  /',  /") . 


(4.9) 


This  solution  differs  from  the  exact  solution  only  in  the  addition  of  the 
damping  correction,  Eq,  (2.  3),  to  the  frequencies. 
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The  vibrations  of  this  nonlinear  system  will  be  considered  as  forced  vibrations 
of  a  linear  conservative  system  with  a  frequency  o>  subjected  to  the  external 
force 


2:f  -  (s*  -  «*)  /  -  2*/'  -  v(/.  /',  n 


Expanding  the  right-hand  term  of  this  expression  in  a  Fourier  series  and 
using  Eq.  (4.  7)  gives 


—  F{a,  to)  sintol  0(a,  to)  costui 


where 


F(a,  w)  —  (to*  —  to*) a  —  0{a,  to), 
Ola,  5)  =  —  2cwa  —  If' (a,  to). 


From  the  beginning  we  have  assumed  that  the  solution  does  not  contain  a 
cosine  part 

Fla,  to)  ==  0. 


This  equation  is  identical  to  the  first  equation  in  Eqs.  (4.  8).  In  the  following, 
we  will  consider  the  equation: 


/"  +  to*/  =--=  fi'ln,  to)  contot. 


(4.  10) 
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If  G(a,  w)  i*  a  slowly  changing  function,  then 


/(I)  =  a(()8ina>( 


where  the  following  equation  is  an  approximate  (asymptotic)  solution  of 
Eq.  (4.  10): 


I 

0 


Differentiation  reduces  the  second  of  Eqs.  (4.8)  to 

In  spite  of  all  their  physical  clarity,  both  justifications  cited  above  are 
not  mathematically  rigorous.  The  rigorous  proof  of  the  method  of  slowly 
changing  amplitudes  was  given  by  L.  I.  Mandel'shtam  and  N.  D.  Papaleksi 
(Ref.  56).  Applying  the  method  of  the  small  parameter,  they  showed  that  the 
first  approximations  correspond  in  both  methods.  The  first  approximation 
according  to  the  method  of  Krylov  and  Bogoliubov  (Ref.  97;  see  also  Ref.  98) 
and  to  the  method  of  Andronov  and  Khackin  (Ref.  99)  lead  also  to  analogous 
results. 

If  higher  harmonics  are  introduced 


/(<)  =  aj(()8ina>(  H-  a,(()sin2(u(  -|-  8in3d>(  +  •  •  •, 


we  can  increase  the  accuracy  of  the  results.  The  smaller  the  nonlinearity 
the  greater  is  the  accuracy  given  by  the  first  harmonic  approximation. 
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4.  The  method  of  slowly  changing  amplitudes  can  also  be  applied  to  the 
problem  of  forced  vibrations.  For  example,  let  us  consider  the  differential 
equation 


/"  +  2*  /'+»•/  +  y(/.  /',  /")  -  S  8in0  f. 


(4. 11) 


We  will  seek  its  solution  in  the  form 

/(<)  3a  a(()>in9  (  +  &(/)  0080  <;  (4.12) 

where  a(t)  and  b(t)  are  slowly  chainging  amplitudes.  Writing  Eq.  (4.  11)  in 
the  form 


r  +  0*/  =  -88in0H-  {0»  -  to*)/  -  2e/'  -  y(/, /'/"),  (4.  13) 

substituting  Eq.  (4.  12)  on  the  right-hand  side  of  Eq.  (4.  13)  and  separating 
out  the  sin  0t  and  cos  0t  terms,  we  obtain 

f"  +  0*/  =  ^’(o,  fr)8in0  ( -f  G{a,  b)  cos0  (  +  -•-.  (4.  i4) 


This  yields 


i’(o,6)  =  d4-(0*-o>*)®-2e06-  0(0,6), 
0(a,6)=  -(0*  -m*)b-2e9a-V{a,b) 


where 


0(o,6)--J-  /  y(/,/'./")»in0l<il. 


in 

f 


V'(a.b)^^f  y(/,/',  ncmBtdt. 


0 
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An  approximate  solution  of  Eq.  (4.  14)  is 


«(0 


I 

0 


I 

''(0  -=  Yg  f  F{r)dt. 

0 


Through  differentiation,  we  obtain  the  van-der-Pol  equations 


dn 

dV 


db 

di 


igFia.b) 


or,  more  explicitly. 


if  =  Y9^- -  2«  Oa  -  '//(«,  6)1, 

Yi-  =  its  +  (0>  -  w»)«  -  2®  06  -  0{rt,  6)]. 


Thus,  instead  of  a  differential  equation  of  second  order,  Eq.  (4.  11), 
we  have  a  simpler  system  of  two  equations  of  first  order  that  are  not  explicitly 
dependent  on  time.  In  the  case  of  stationary  vibrations 


For  determining  the  steady-state  amplitudes,  we  obtain  the  following  system 
of  algebraic  equations: 


--  (0*  -•  0}*)  b  —  2e  9  a  —  '/'(« ,  6)  =-  0, 
5  ;■  (O’  —  a>*)  a  —  2e  0  b  —  0(rt,  6)  —  0. 


(4.  15) 
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This  equation  system  can  also  be  derived  by  other  means.  For  example,  one 
can  substitute  Eq.  (4.  12)  directly  into  Eq.  (4.  11),  which  will  be  written  in 
the  form  L(f,f',f')  =  0.  One  then  requires  that  the  result  of  the  substitution 
(in  accordance  with  the  variational  method  of  Galerkin)  be  orthogonal  to  the 
"functional  coordinates"  sin  $t  and  cos  9t: 


•  t 

f  L(f,  n  anOtdt  =  0,  /  Lit,  /',  n  coaetdi  -  0. 

0  0 

This  method  directly  leads  to  Eqs.  (4.  15). 

•17.  FREE  VIBRATIONS  OF  A  NONLINEAR  SYSTEM 

1.  Let  us  consider  the  problem  of  the  free  vibrations -of  a  rod,  taking  into 

2 

account  the  nonlinear  terms  introduced  in  the  previous  chapter.  The 
differential  equation  of  this  problem  is 


/"  +  2c/'4-w*/  +  V(/./'.n  =  0: 


(4.  16) 


where  4'(fif'if")  is  a  nonlinear  function  of  displacements,  velocities  and 
accelerations.  Limiting  ourselves  to  quantities  of  less  than  third  order,  the 
nonlinear  function  can  be  represented  in  the  form 

vifj'.n  =  yn  +  +  2xi[fr  +  im  (4.  i?) 


The  first  term  considers  the  influence  of  nonlinear  terms  of  static 
origin,  the  second  term  considers  the  nonlinear  character  of  damping,  and 


A  comprehensive  theoretical  and  experimental  investigation  of  the  free 
vibrations  and  forced  vibrations  of  a  nonlinear  system  can  be  found  in 
Refs.  64,  65,  66,  and  67. 
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the  third  term  considers  the  influence  of  the  inertia  forces  that  arise  from 
longitudinal  displacements.  In  the  following  paragraphs  we  will  discuss  non¬ 
linear  elasticity,  nonlinear  damping,  and  nonlinear  inertia. 

It  is  easy  to  show  that  the  differential  equation  of  the  problem  belongs 
to  a  class  of  differential  equations  with  small  nonlinearities.  This  charac¬ 
teristic  follows  from  the  problem  itself.  As  is  well-known,  the  linear  treat¬ 
ment  of  problems  on  the  free  vibrations  of  rods  gives  results  that  are  in 
good  agreement  with  experimental  data.  The  additional  nonlinear  terms  in 
Eq.  (4.  16)  are  a  correction  through  which  the  linear  approximation  is  made 
more  precise. 

To  obtain  a  quantitative  estimate,  one  must  reduce  Eq.  (4.  16)  to  a 
nondimensional  form  and  apply  the  condition  in  Eq.  (4.  2).  The  calculation 
will  be  left  to  the  reader.  During  stable  equilibrium  the  influence  of  non- 
linearities  depends  upon  the  amplitude  of  the  vibrations  fits  growth  is  pro¬ 
portional  to  the  square  of  the  amplitude).  This  means  that  the  vibration 
amplitudes  must  be  sufficiently  small.  In  all  practical  problems  this  condi¬ 
tion  is  assumed  to  be  satisfied. 

2.  We  will  now  consider  calculations.  Substituting  the  expression 
f  =  a  sinwt  into  Eq.  (4.  17)  and  introducing 

ain'col  —  sin3ci><)> 

sina>lo(M*wl  =  i(ain«»<  -f  MinSoil), 

8in*wlcoad>l  =,  -^■(oo«a>l  —  oosSoif) 


we  obtain 


y(/i  /',  /")  =  —  N9ja*iinul  -f  -^-a*oo»at  +  ■  •  • . 


(4. 18) 
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In  accordance  with  the  method  of  slowly  changing  amplitudes,  the  derivatives 
a'  and  a",  and  also  terms  containing  higher  harmonics,  are  neglected.  From 
Eq.  (4.  18)  it  follows  that 


0(a,  co)  K  _  Kai*ja*, 

so  that  Eq.  (4.  8)  for  our  case  takes  the  form 


and 


«»•  —  +  —  xc5*Jo*™0 


(4.  19) 


(4.  20) 


Equation  (4.  19)  permits  one  to  determine  the  frequency  of  the  free 
vibrations  of  the  nonlinear  system 


a>  =  (V 


f 


4  w* 

1  +  v«*  ■ 


(4.21) 


As  one  can  see  from  Eq.  (4.  19),  the  natural  frequency  of  the  nonlinear  system 
depends  upon  the  vibration  amplitude.  Therefore,  the  existence  of  nonlinear 
elasticity  leads  to  an  increase  of  the  frequency  with  amplitude;  conversely, 
nonlinear  inertia  causes  a  decrease  of  the  natural  frequency.  If  we  consider 
that  the  nonlinearity  is  small,  we  can  write  Eq.  (4.21)  in  the  form 
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For  X  <  {3/4)('y/w^)  the  frequency  of  free  vibrations  will  increase  with 
a  growth  of  amplitude;  for  x  >  (3/4)('y/w^)  the  frequency  will  decrease  with 
an  increase  of  amplitude.  For 


(4.23) 

the  effect  of  nonlinear  elasticity  and  nonlinear  inertia  is  the  same;  the  free 

vibrations  remain  isochronous.  In  this  sense  a  nonlinear  inertia  with  a  co- 

efficient  X  and  a  nonlinear  elasticity  with  a  coefficient  y  =  (4/3)X<*>  are 

equivalent.  We  are  reminded  that  an  analogous  correspondence  arises  for 

the  linear  system:  the  influence  of  transverse  forces  of  inertia  is  equivalent 

to  the  influence  of  a  continuous  elastic  foundation  with  a  resistance  coefficient 
2 

of  k  =  -mw  .  The  factor  4/3  in  our  case  can  be  explained  by  the  nonlinear 
character  of  the  system. 

Let  us  consider  the  case  when  the  longitudinal  elastic  connection  is 
absent.  Then,  as  we  will  show,  the  nonlinear  inertia  term  always  turns  out 
to  be  dominant  in  this  case.  The  coefficient  of  nonlinear  elasticity  according 
to  Eq.  (3.  16)  is 


y= 


8P  • 


(4.  24) 


Since  the  imperfect  elastic  material  usually  gives  a  "soft"  nonlinearity,  one 
must  consider  Eq.  (4.  24)  as  the  maximum  value  of  the  coefficient  of  nonlinear 
elasticity  in  the  absence  of  nonlinear  connections.  On  the  other  hand,  the 
smallest  value  of  the  coefficient  X  given  in  Eq.  (3.  20)  is 


X  = 


41*  \3 
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{ 


( 


If  we  consider  that  for  this  value 


4  NIB* 

3  y 


ts*  /  1  3  \ 

=  “jriy“87?; 


m  8 


we  arrive  at  a  conclusion  which  is  irhportant  for  all  future  considerations. 

In  the  absence  of  longitudinal  couplings  the  predominant  nonlinear  factor  is 

nonlinear  inertia.  This  conclusion  is  valid  for  all  rods  which  con^fetitute 

statically  determinate  systems  and  also  to  those  rods  of  statically  indeter- 

2 

minate  systems,  for  which  X  >  (3/4)(y/(»)  ). 

Let  us  now  estimate  the  order  of  magnitude  of  the  nonlinear  correction 
to  the  frequency.  If  we  substitute  the  values  of  coefficients  y  X  into 
Eq.  (4.  22),  we  can  write  this  equation  in  the  form 


where 


k 


S.2  +  1S-^J. 


denotes  the  "longitudinal  mass"  which  is  fixed  to  the  moving  end  of  the 
rod.  The  mass  of  the  rod  is  not  considered  during  computations.  (It  is 
accounted  for  in  the  term  3.  2. ) 

Let  M.  =  0.  In  this  case  the  change  in  the  free  vibration  does  not  ex- 
ceed  3  percent  even  at  such  a  large  amplitude  as  a/i  =  0.  1.  If  the  rod  is 
part  of  a  rod  system,  the  influence  of  the  nonlinear  terms  can  then  increase 
by  ten  or  one  hundred  times.  This  can  be  shown  by  examining  the  example 
in  sl4. 
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An  oscillogram  obtained  experimentally  of 
the  vibrations  of  a  rod  with  a  large  inertia  non¬ 
linearity  is  shown  in  Fig.  23.  The  change  of 
period  of  the  free  vibration  can  easily  be  seen. 

3.  For  an  explanation  of  the  law  of  damping  one 
must  return  to  Eq.  (4.  20).  If  we  write  it  in  the 
form 


edi  = 


.d» 


we  obtain  upon  integration 


—el  =  Ina  —  y 


or 


e-'*. 


The  constant  integration  is 


w 

c  =  l-^ 


+  <4 


where  a^  is  the  initial  amplitude.  Therefore 
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Solving  this  equation  with  respect  to  the  amplitude  a,  we  finally  obtain 


a  =5 


tt 

4( 


(4.  25) 


For  the  special  case  of  linear  damping, 


«  =  (4.26) 

Nonlinear  damping  gives  a  more  rapid  decrease  of  the  amplitudes,  as  one 
might  expect. 

For  practical  purposes  the  following  interpretation  of  Eq.  (4.  20)  might 
be  useful.  If  one  divides  the  whole  damping  period  into  sufficiently  small 
intervals  of  time,  one  can  average  the  nonlinear  part  of  the  damping  by 
assuming  it  is  constant  during  every  interval; 


(4.27) 


Here  a  is  an  average  value  of  the  amplitude  in  the  given  time  interval.  For 
every  interval  Eq.  (4.  20)  takes  the  form 


da 

-di  = 


from  which 


a  =  0(6“** 
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Thus  Eq.  (4.  20)  describes  damped  vibrations  that  occur  with  a  variable 
damping  constant.  Such  a  "quasi -linear"  treatment  is  all  the  more  suitable 
since  essentially  it  gives  the  basis  for  a  generally  accepted  method  of  re¬ 
ducing  experimental  diagrams  of  damped  vibrations.  As  is  well-known,  the 
experimental  damping  constant  is  determined  from  the  condition 

‘  =  (4.28) 

whether  it  is  constant  during  the  entire  process  of  damping  or  changes  with 
amplitude.  In  the  last  case  the  damping  is  nonlinear,  or  more  precisely,  it 
does  not  follow  the  exponential  law  of  Eq.  (4.  26).  The  application  of  Eq.  (4.  28) 
assumes  that  the  damping  is  "linearized"  within  the  "limits"  of  each  interval 
At. 


#18.  FORCED  VIBRATIONS  OF  A  NONLINEAR  SYSTEM 

1.  We  will  discuss  briefly  the  problem  of  forced  vibrations  of  a  nonlinear 
system.  ^  The  steady-state  amplitudes  of  these  vibrations  can  be  determined 
from  Eq.  (4.  15).  The  amplitude  of  the  "generalized  force"  is  denoted  by  S, 
which  will  be 


5  = 


2P 


for  example,  in  the  case  of  a  rod  loaded  at  the  center  (Fig.  24). 


Fig.  24 


Substituting  the  approximate  solution 


f{t)  - a  sill  9  (  —  b  cos  0  t 


See  footnote  2  on  page  41. 
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( 


into  the  expression  for  the  nonlinear  function,  Eq.  (4.  17),  and  performing 
the  operations  indicated  in  the  preceding  paragraph,  we  obtain 


4Ha.b)^A*  (Ay  _k0*)«  +  ^06]. 
!P(a,  fc)  =  -  (Ay  _  X  »•)  6  + -^ Ott]. 


2  2  2 

Here  A  is  the  amplitude  of  the  vibrations,  i.  e.  ,  A  =  a  +  b  .  Therefore 
Eqs.  (4.  15)  take  the  form 


—  (9*  —  «»*)6  —  2e9«  —  A*|-(Ay  —  =  0. 

«  +  (flt  _  o,t),  _  2,96  _^i[(Ay  _  *e.)«  +  ^96)  =  0 


(4.  29) 


The  equation  system,  Eq,  (4.  29)  is  unsolvable  in  the  general  form; 
therefore  we  restrict  ourselves  to  the  case  in  which  the  damping,  either 
linear  or  nonlinear,  can  be  neglected.  The  equation  system  for  this  case  is 
simplified  to 


-(0* -  «•) 6  +  (Ay  -  s fl*)  =  0, 
8  +  {9*—  o>*)«  —  A*a  ( ^  y  —  ® 


and  obviously  can  be  satisfied  by  b  =  0,  A  =  a.  The  amplitude  of  the  steady- 
state  vibrations  is  determined  from  the  cubic  equation 


«  +  (9»  -  w*)  d  -  i4*(Yy  -  x9*)  =  0. 


(4.  30) 
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2.  The  roots  of  Eq.  (4.  30)  can  be  graphically  determined  as  the  coordinates 
of  intersection  points  of  a  straight  line 


(4.31) 


with  the  cubic  parabola  y  =  A  .  Figure  25  corresponds  to  the  case  in  which 

2 

the  nonlinear  elasticity  is  dominant,  i.  e.  ,  (3/4)y  -  Xu  ^  0- 


During  changes  of  the  exciting  frequencies  from  zero  to  infinity,  the 
straight  line,  Eq.  (4.  31),  rotates  clockwise  from  the  initial  position  MN  to 
the  end  position  M^N^.  As  long  as  the  frequency  is  sufficiently  small  (at 
least  for  0  <  u),  the  straight  line  intersects  the  cubic  parabola  once,  i.  e.  , 

Eq.  (4.  30)  has  one  real  root.  For  larger  frequencies,  three  real  roots 
occur  and,  consequently,  three  possible  amplitude  values  also  occur.  At 

(4.32) 
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the  straight  line  takes  a  horizontal  position  so  that  Eq.  (4.  30),  in  addition  to 
having  the  one  finite  root,  has  also  the  roots  A  =  ±oo. 


The  resonance  curve  is  shown inFig.  26. 

Its  character  is  different  not  only  from  the 

resonance  curve  of  the  linear  problem  but  also 

from  the  resonance  curve  for  the  case  in  which 

4 

only  nonlinear  elasticity  is  present.  As  in 
the  much  investigated  case  of  nonlinear  elas¬ 
ticity,  from  the  three  correct  solutions  which 
exist  in  the  region  CD,  only  two  will  be  stable. 
The  unstable,  i.e.  ,  the  physically  impossible 
solution,  is  denoted  by  a  dotted  line. 


The  existence  of  two  stable  solutions  ^6 

leads  to  a  result  which  is  called  the  "over¬ 
hang."  By  a  gradual  increase  in  the  excitation  frequency  we  can  bring  the 
system  into  the  region  where  two  stable  solutions  exist  (Fig.  27a).  The 
growth  of  the  amplitudes  will 
proceed  initially  along  the  curve 
until  vibrations  occur  at  a 
certain  point  M  j .  The  amplitude 
(not  calculated  with  respect  to 
time)  decreases  suddenly  to  the 
magnitude  MM2f  and  then  de¬ 
creases  along  the  curve 
By  reversing  this  frequency 
process,  the  amplitude  of  the 
.  vibrations  will  initially  increase 
smoothly  along  the  curve  l-'2^2 
(Fig.  27b).  The  sudden  change 


This  last  case  examined  in  many  textbooks  on  the  theory  of  vibrations 
(for  example.  Ref.  98)  can  be  obtained  from  our  equations,  if  in  them  one 
sets  X  =  0.  The  values  of  the  amplitude  increase  for  increasing  frequencies. 
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now  occurs  at  point  N2.  after  which  the  amplitudes  follow  the  curve  NjKj. 
This  shows  that  the  amplitude  of  the  forced  vibrations  depends  not  only  on 
excitation  frequency  in  the  considered  instant  of  time  but  also  on  the  "history" 
of  the  nonlinear  system. 


The  question  of  the  "overhang"  will 
be  examined  in  connection  with  parametrically 
excited  vibrations  (eZB). 

We  have  examined  the  case  in  which 
the  nonlinearities  are  due  essentially  to  elas¬ 
ticity.  If  the  influence  of  nonlinear  inertia  is 
predominant,  the  resonance  curve  bends 
toward  the  side  of  decreasing  frequencies; 
the  "overhang"  of  the  resonance  curve  also 
bends  towards  the  side  of  decreasing 
frequencies  (Fig.  28), 
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CHAPTER  FIVE 


AMPLITUDES  OF  VIBRATIONS  AT  THE  PRINCIPAL 
PARAMETRIC  RESONANCE 


•19.  BASIC  EQUATIONS 
1.  Let  us  investigate  the  differential  equation 


r  +  2f  /'  -f-  l?»(l  -  2/4  cm®  0  /  t-  fil,  /'. !")  =  0. 


(5.  1) 


that  describes  parametrically  excited  vibrations  including  the  effects  of  non¬ 
linear  factors,  and  let  us  formulate  the  problem  of  finding  the  solutions 
that  correspond  to  steady-state  vibrations. 

It  was  established  in  *8  that  the  solutions  of  the  linear  problem  have 
the  form 


(5.2) 


on  the  boundaries  of  the  first,  third,  and  generally,  the  odd  regions  of  dynamic 
instability  where  a^^  and  b^^  are  constant  coefficients.  It  would  be  natural  to 
seek  the  solution  of  the  nonlinear  problem  in  such  a  form  that  is  valid  within 
the  limits  of  the  odd  regions  of  instability.  The  following  considerations  lead 
to  this  form  of  solution. 


For  a  specific  choice  of  its  coefficients,  Eq.  (5.  2)  can  satisfy  Eq.  (5.  1). 
In  fact,  the  result  of  the  substitution  of  the  series  into  Eq.  (5.  1)  will  not  con¬ 
tain  other  periodic  terms  except  sin  (k9t)/2  and  cos  (k9t)/2  for  odd  values  of 
k.  However,  this  applies  only  if  the  nonlinear  function  i|4(f ,  f ,  f ”)  does  not 
contain  terms  of  even  powers.  The  function 

y  (/ .  /")  -  y  /•  +  2t,.  /*  /'  +  2ic  /[/  r  +  (/')*]  (5.3) 

will  satisfy  this  requirement. 
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The  third  and  deciding  reason  for  the  use  of  Eq.  (5.  2)  is  the  experi¬ 
mental  results,  which  indicate  that  the  steady- state  vibrations  within  the 
limits  of  the  odd  regions  of  instability  have  exactly  such  a  form.  We  will 
postpone  a  discussion  of  the  experimental  data  until  s22. 

2.  Let  us  seek  the  solution  of  Eq.  (5.  1)  in  the  form  of  Eq.  (5.  2).  Substi¬ 
tution  gives 


-  +  2)  ^  +  00.(4  -  2)  Y]  +  y*(/.  n  =  0. 


The  nonlinear  function  in  Eq.  (5.  4)  is  expanded  in  a  Fourier  series 


¥Uf.  r.  n  -  1'  {0tn\n^L  +  Vtcon^\ 

^  I*  3*  5  '  * 


2  r 


(5.  5) 


Its  coefficients 


4m 

•  «* . y*(/.  /'.  n "in  *J‘  dl. 

0 

4j* 

. ^  /  y*  (/,  r.  n  «>H  *•-  dl 

0 


(5.6) 


are  obviously  nonlinear  functions  of  the  coefficients  of  Eq.  (5.  2)  and  in  the 
case  of  Eq.  (5.  3)  are  homogeneous  functions  of  third  order.  Substituting 
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Eq.  (5.  5)  into  Eq.  (5.4)  and  equating  the  coefficients  of  sin  (k9t)/2  and 
cos  (k9t)/2,  we  obtain  the  following  system  of  equations 


+n)—  **|«i  —  rOb,  —  fci)  =  0, 

[/3»(1  -  /I)  -  ^  ]  ft,  +  e  0«,  -  +  V,(ai .  6,)  =  0, 

—  teOht  —  +  «»*i)  +  ^()  =  0> 

^ -f  ktOa^  —  /7*/i(6t-t  +  bk-n)  +  ft()  =  0, 

(fc  =  3.8,...). 


(5.7) 


For  brevity  we  write 


•  •  ••  fc|.  ^»»  •  •  •)  =  *()• 

V'iiat ft| .  fc».  • .  •)  =  *'*(*1.  ^<). 

3.  In  order  to  determine  the  amplitudes  of  the  steady-state  vibrations 
within  the  limits  of  the  even-numbered  regions  of  excitation,  we  will  seek 
solutions  in  the  form 


/(<)  =  6, ^  ^ sin  +  ftt  .  (5.8) 

This  form  is  suggested  by  considerations  analogous  to  those  stated  in  No.  1. 

Substituting  Eq.  (5.  8)  into  Eq.  (5.  1)  and  expanding  the  nonlinear  func¬ 
tion  in  the  series 


y*{/.  r,  n  =  y',  + 

* 


V 
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we  obtain  the  following  system  of  equations; 


{Q*  -  0*)rt,  -  -  2rf;fcj  -I  h,)  r.  o, 

(fj*  -  0*)fc,  -  i»!J*V2h,  !  fc,)  I  -.‘f  Off,  -1  y^,K.  ft,) .  0. 

j«»  -  /«//*(«*-•  4  n*  ;•)  —  if  Oft,  i-  «/»,(»!(,  fti)  =-  t», 

(O*  -  i  +  ft,,»)  i  i  f  Off*  •//*(ff,,  6.)  .  0 

(i=  4, 


(5.9) 


Coefficients  and  ^  for  k>2  are  calculated  according  to  Eq.  (5.6),  and 
the  coefficient  according  to  the  formula 


4m 

«l . .  ■  ■)  -  4®,  f  V*U,  /'•  /")rf< 

A 

We  note  in  conclusion,  that  Eqs.  (5.  7)  and  (5.  9)  could  be  obtained 
differently,  i.e.  ,  by  means  of  the  Galerkin  variational  method.  If  one  sub> 
stitutes,  for  example,  Eq.  (5.2)  into  Eq.  (5.  1),  one  can  require  that  the 
corresponding  expression  of  Li(f,  f',f  ")  be  orthogonal  to  each  of  the  functional 
coordinates  sin  (k9t)/2  and  cos  (k9t)/Z: 


f  t 

;  **  " 
(i  =. 


These  expressions  reduce  to  Eq.  (5.  7). 
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•  20.  DETERMINATION  OF  STEADY -STATE  AMPLITUDES* 

1.  If  we  investigate  the  vibrations  for  the  principal  resonance  6  ZH, 
we  can  neglect  the  influence  of  higher  harmonics  in  Eq.  (5.  2)  and  can  assume 


/(f)  =  a  Mill  h  c«>«-  -  - 

z  z 


(5.  10) 


as  an  approximation.  Therefore,  Eqs.  (5.  7)  are  essentially  simplified, 
and  a  system  of  two  equations  for  the  coefficients  a  and  b  remains: 


|/A(I  f  /d  -  *’|a  -  -t-  </*(«,  h)  ...  0, 

j/A(l  --  ft)  -  j*  1  rfln  -i-  '/'(«,/>)  -  0. 


(5.  11) 


For  the  determination  of  the  quantities  ^a,b)  and  ^'(a,b),  we  sub¬ 
stitute  Eq.  (5.  10)  into  Eq.  (5.  3).  Computations  give 

V* (/•/'. /")  -  r  1,0 h  - 

*■  ®  I  *1.  ® “  ^0*  h)  iKM  !■•••. 


Terms  containing  higher  harmonics  are  not  written  out;  the  amplitude  of 
steady- state  vibrations  is  denoted  by  A  where 


Thus,  the  first  coefficients  of  the  expansion  of  the  function  i;i*(f,  f,  f")  in 
a  Fourier  series  are 

Jsee  Ref.  57. 

The  method  of  the  slowly  changing  amplitude  also  leads  to  this  approximation. 
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-  ft®* -»«•*•). 

¥>(,.  *)  =.  if  (Syfc  +  «*®a  -  K0*b). 


(5.12) 


2.  For  further  investigation,  we  will  write  the  system  of  Eqs.  (5. 11)  in 
the  following  form: 

(1  +  -  ..••)  -  0,  j 

(I  _  ^  _  „i)6  +  +  .!*(  fc  +  "f*-®  -  MS*b)  «  0 

where 


J  •=  -  rt*- .  '■ 


(5.14) 


It  is  obvious  that  Eqs.  (5.  13)  will  be  satisfied  for  a  =  b  =  A  =  0.  This 
solution  corresponds  to  the  case  when  transverse  vibrations  of  the  rod  are 
absent. 

One  can  find  the  non>zero  solutions  in  the  following  manner.  We  will 
consider  Eq.  (5.  13)  as  a  system  of  homogeneous  linear  equations  with  respect 
to  a  and  b.  This  system  has  solutions  that  differ  from  zero  only  in  the  case 
when  the  coefficients  of  the  determinant  disappear; 

1  +  #»  -  s*  —  ii*  I* a*  —  “  (d  +  dtd») 

I  (A  +  At  A*)  I  -  ^  -  a»  -  d*  (a  »•  - 

Expanding  the  determinant  and  solving  the  equation  obtained  with  respect  to 
the  amplitude  A  of  the  steady- state  vibrations,  we  find 


=  0.  (5.15) 


A 


r 


,(1  _  ,•)  (j»*+  ^l)  -  "I  fj*d  +  did  -  a»)l* 


(5.16) 
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where 


p  •=  HU* 


8y 

riF- 


(5.  17) 


3.  Equation  (5.  16)  is  too  cumbersome  for  further  investigation.  Therefore 
we  shall  first  clarify  how  the  character  of  the  resonance  curves  depends 
upon  the  form  of  the  nonlinear  function.  For  the  case  t)'  case  of 

nonlinear  elasticity),  Eq.  (5.  16)  gives 

A  diagram  of  the  dependence  of  the  amplitude  on  the  frequency  is 
shown  in  Fig.  29a,  where  the  two  plotted  solutions  correspond  to  the  two 
signs  in  Eq.  (5.  18).  One  of  these  solutions 
(represented  by  the  dotted  line)  is  obviously 
\mstable.  It  is  characteristic  for  the  case  of 
nonlinear  elasticity  that  resonance  curves  are 
bent  towards  the  side  of  the  increasing  exciting 
frequencies.  We  have  already  encountered 
analogous  properties  during  the  investigation 
of  forced  vibrations  (  #18). 

Let  us  consider  the  case  of  nonlinear  damping 


a) 


b) 

Fig.  29 


c) 


Equation  (5.  16)  assumes  the  form 


A  - 


nA 

n 
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The  corresponding  resonance  curve  is  represented  in  Fig.  29b.  As  is  seen 
from  the  figure,  the  resonance  curve  in  this  case  is  approximately  sym¬ 
metrical  with  respect  to  the  maximum. 

Finally,  for  the  case  of  the  nonlinear  inertia 

y(/. 


we  obtain  the  amplitude  equation 


(5.  19) 


Here  we  also  have  two  solutions,  one  of  which  is  unstable  (Fig.  29c).  Con¬ 
trary  to  the  case  of  nonlinear  elasticity,  the  resonance  curves  are  bent 
toward  the  side  of  decreasing  exciting  frequencies. 

We  note  that  the  base  of  the  resonance  curves  in  all  three  cases  (as  in 
the  general  case)  does  not  depend  on  the  magnitude  and  the  character  of  the 
nonlinearity  but  coincides  with  the  interval  of  the  instability  determined  by 
the  methods  of  linear  theory.  In  fact,  if  we  set  A  =  0  in  Eq.  (5.  15),  we 
obtain 
% 

+  -if  I  . 


coinciding  with  the  equation  of  the  critical  frequencies,  Eq.  (2.  12). 

It  is  evident  that  in  the  case  of  nonlinear  elasticity  or  nonlinear  inertia, 
it  is  possible  to  have  an  "overhang"  of  the  resonance  curve  beyond  the  limits 
of  the  region  of  excitation.  (Compare  with  *18.  )  In  the  first  case,  the 
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"overhang"  will  occur  towards  the  side  of  increasing  frequencies;  in  the 
second  case,  it  will  occur  towards  the  side  of  decreasing  frequencies.  This 
phenomenon  will  be  investigated  in  detail  in  *23. 

4.  The  form  of  the  resonance  curves  for  each  of  the  three  fundamental 
cases  is  so  typical,  that  according  to  the  character  of  the  curves,  it  is  pos¬ 
sible  to  judge  which  nonlinear  factors  predominate  in  each  case.  It  can  be 
seen  from  Eq.  (5.  16),  that  the  manner  in  which  the  resonance  curve  is  bent 
depends  on  the  sign  of  the  quantity 


For 


4  »tfl« 

a  r 


>1 


(»•!) 


the  resonance  curves  bend  toward  the  side  of  decreasing  frequencies,  i.  e.  , 
the  amplitudes  will  increase  as  the  exciting  frequencies  decrease.  If, 
howe  ve  r , 


<1 


(s  <•  1), 


the  resonance  curves  will  then  bend  toward  the  side  of  increasing  frequencies. 

As  in  the  case  of  free  vibrations  (  *17),  it  can  be  shown  that  if  the  rod 
does  not  have  longitudinal  elastic  springs,  then  the  nonlinear  inertia  is  the 
decisive  nonlinear  factor.  From  the  relation  derived  in  *17 

3  ~  y  *^3  ' '3  8»M 
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it  follows  that 


4  nn^O* 

T~T~ 


Sj.* 

3 


f "  \JL 

U  iJTj  4«i*  • 


However,  for  the  principal  region  of  instability 


S«*  /  I 

"riT 


> 


1, 


only  if  the  longitudinal  force  is  not  in  the  vicinity  of  the  Euler  value. 

The  dominating  influence  of  nonlinear  inertia  is  confirmed  by  experi¬ 
mental  data.  Oscillograms  of  parametric  excited  vibrations  taken  near  the 
principal  resonance  $  =  ZQ  are  shown  in  Fig.  30.  Thede  oscillograms  show 


Fig.  30 


the  vibrations  which  arise  for  a  decreasing  exciting  frequency.  The  first 
frame  corresponds  to  the  upper  boundary  of  the  region  of  parametric  excita¬ 
tion  (0  =  42.  0  cps).  Figure  11a  is  an  enlargement  of  this  frame.  Along  with 
the  parametric  excited  vibrations,  forced  vibrations  also  occur  with  a  doubled 
frequency  (i.e.  ,  with  a  frequency  of  the  external  force).  The  additional  frames 
correspond  to  a  gradual  decrease  of  the  exciting  frequency.  It  is  charac¬ 
teristic  that  the  accompanying  vibrational  amplitude  increases.  The  sixth 
frame  corresponds  to  an  "overhang"  of  the  resonance  curve  (0  =  36.  5  cps), 
after  which  the  "break"  of  the  amplitude  occurs.  Only  forced  vibrations 
having  a  frequency  of  the  external  force  remain. 
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The  dependence  of  the  amplitude  on  the 
frequency  also  can  be  seen  in  Fig.  31.  In  con¬ 
trast  to  the  oscillogram  in  Fig.  30,  this  oscil¬ 
logram  shows  a  continuous  record.  Adecrease 
of  the  exciting  frequency  by  20  percent,  which 
takes  the  system  beyond  the  limits  of  the  region 
of  instability  along  the  "overhang,  "  increases 
the  amplitudes  by  a  factor  greater  than  3. 

e21.  INVESTIGATION  OF  THE  EQUATION  FOR 
STEADY -STATE  AMPLITUDES 

1.  We  will  first  investigate  the  influence  of  linear  damping  on  the  magni¬ 
tude  of  the  steady-state  amplitudes.  For  this  purpose,  we  neglect  in 
Eq,  (5.  16)  the  terms  which  take  into  account  the  nonlinear  portion  of  the 
damping.  This  can  be  represented  in  the  following  form 


summ  "H  h 


(5.  20) 


where  as  before 


Equation  (5.  20)  gives  real  values  for  the  amplitude  in  the  case  (nA)/‘ir<  p,  or 
because  n»l,  in  the  case  d/ir<|ji.  If  however,  then  steady- state 

vibrations  do  not  arise.  This  result  is  found  to  be  in  complete  agreement 
with  the  results  of  linear  theory. 

On  the  other  hand,  for  sufficiently  small  damping,  its  influence  on  the 
magnitude  of  the  amplitudes  is  imperceptible  for  practical  purposes.  For 
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example,  at  n  =  1,  Eq.  (5.20)  assumes  the  form 


A 


where  is  the  critical  value  of  the  excitation  parameter  (#9).  Even  at 
Pi)t  =  l/2p,  the  influence  of  damping  comprises  an  order  of  magnitude  of 
6  percent.  With  a  further  decrease  of  damping,  the  approximate  values 
which  were  computed  for  the  conservative  case  are  still  more  exact.  This 
is  seen,  for  example,  from  Fig.  32,  which  represents  the  dependence  of  the 

dimensionless  amplitude  on  the  ratio  of  the  excitation 
parameter.  In  general  at 


>*>®/**  (5.21) 


Fig.  32  the  simple  formula  of  the  conservative  problem 


(5.  22) 


gives  sufficient  accuracy  for  practical  purposes.  In  particular,  the  largest 
amplitudes  within  the  limits  of  the  region  of  instability  (i.  e.  ,  without  the 
consideration  of  the  "delay")  can  be  determined  according  to  the  formula 


A  ~ 


(5.  23) 


In  the  case  of  predominant  nonlinear  inertia  (p  >  0),  the  greatest 
amplitudes  are  attained  on  the  lower  boundary  of  the  resonance  region;  in 
the  case  of  predominant  nonlinear  elasticity  (p  <  0),  the  greatest  amplitudes 
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ftre  attained  on  its  upper  boundaries.  Equation  (5.  22)  can  also  be  written  in 
the  following  form; 


A* 


(5.  24) 


Here,  the  upper  or  lower  critical  frequency  is  taken,  depending  on  the  sign 
of  p  instead  of  For  the  proof,  let  us  investigate  Eq.  (5.  15)  at 

A  =  =  0: 


l+ft-n*-A*p  0 

0  I  -  p  -  m*- A*  p 


0. 


If  we  compare  it  with  Eq.  (2.  12)  for  the  frequency  boundaries,  we  see  that 
2  2  2 

it  is  satisfied  at  n  -t-  A  p  =  n^  .  Consequently, 


follows  from  Eq.  (5.  24). 

In  what  follows,  unless  otherwise  stipulated,  we  will  assume  that  the 
condition  in  Eq.  (5.  21)  is  fulfilled,  and  we  will  make  use  of  Eqs.  (5.  22), 

(5.  23),  and  (5.  24). 

2.  We  shall  now  turn  to  the  question  concerning  the  influence  of  the 
magnitude  of  external  loading  on  the  steady- state  amplitudes. 

As  seen  from  Eq.  (5.23),  the  amplitude  of  the  vibrations  grows  in 
proportion  to  the  square  root  of  the  excitation  parameter 

^  =  (5.25) 
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i.  e.  ,  proportional  to  the  square  root  of  the  amplitude  of  the  periodic  force. 

A  nonlinear  dependence  between  the  loading  and  the  vibration  amplitude  is 
usually  characteristic  for  problems  which  are  described  by  nonlinear  differ¬ 
ential  equations,  However,  the  amplitude  of  the  longitudinal  displacements 
at  the  point  of  application  of  the  force  also  grows  in  proportion  to  the  excitation 
parameter. 

At  first  glance,  the  constant  component  of  the  longitudinal  force 
must  increase  the  amplitudes  of  the  vibrations.  Actually  an  increase  of 
the  force  also  increases  the  excitation  parameters: 


A 

HP*  -  Pt)  ■ 


However,  this  dependence  is  more  complicated.  The  force  Pq  is 
usually  of  gravitational  origin,  i.  e.  ,  one  way  or  another  it  is  related  to  the 
weight,  and  therefore  an  increase  usually  causes  an  increase  in  the  non¬ 
linear  inertia  of  the  system.  This  also  can  bring  about  a  decrease  of  the 
amplitudes.  We  will  illustrate  thii  by  a  simple  example  (Fig.  20). 

Let  the  simply- supported  rod  be  loaded  by  the  force  Pq  -f  P^  cos  0t, 
where  the  force  Pq  is  associated  with  the  rod  mass  ^Q/g*  Neglecting  the 
influence  of  nonlinear  elasticity  and  linear  and  nonlinear  damping,  we  obtain 

^  n*{P,+  ka). 

401* 

where  G  =  mgi  is  the  free  weight  of  the  rod,  and  k  is  the  coefficient  in 
Eq.  (3.20)  designated  by  X-  Fq.  (5.22)  at  d  =  2R  gives  A^  =p/X,  or 

Jl»  AG 

^  ~  {P,  +  tG)iPt-  P^  '  (5.26) 
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The  dependence  of  the  amplitude  on  the  constant  component 
Pq  is  shown  in  Fig.  33.  At  Pq  =  0.  the  vibration  amplitudes 
are  relatively  large: 


ifp, 

“  »•*  /*, 


An  increase  of  the  force  Pq  decreases  the  amplitudes  until 
a  minimum  is  reached  at  Pq  =  1/2(P  H-  kG)  »  (1/2)P^.  Subsequent  amplitudes 
once  again  increase. 


3.  In  the  presence  of  two  factors,  i.  e.  , 
nonlinear  inertia,  the  resonance  curves  as 
Figure  34a  corresponds  to  the  case  p  >  0 
(predominant  nonlinear  inertia). 

As  seen  from  the  diagrams,  the 
"overhang"  in  this  case  is  bounded  by  the 
frequency  0^.  This  frequency  is  deter¬ 
mined  by  equating  the  expression  for  p 
to  zero: 


XN' 


_ 

40* 


=  0. 


the  nonlinear  elasticity  and  the 
sume  the  form  shown  in  Fig.  34. 


Hence, 


(5.  27) 


The  magnitude  ^  y/x  ie  none  other  than  the  linearized  frequency 

of  the  "longitudinal  system"  (the  nonlinear  elasticity  plus  the  nonlinear 
inertia).  Thus  in  the  case  of  an  end  mass  M.  and  a  connection  with  a  spring 
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constant,  we  have 


k  = 


4m  P 


r.  y  = 


it*e 
4m  P 


which  gives  u. 


c/M.  .  Equation  (5.  27)  can  be  interpreted  as  the  condition 
Li 


for  the  occurrence  of  resonance  in  the  "longitudinal"  system  although  it 
differs  from  the  condition  of  synchronization  by  the  factor  of  for  the 
natural  frequency. 

If  the  nonlinear  elasticity  and  the  nonlinear  inertia  compensate  one 
another,  then  the  frequency  9^  lies  in  the  region  of  dynamic  instability. 

This  is  the  most  unfavorable  case  from  the  point  of  view  of  the  amplitudes 
of  the  vibrations:  the  parametric  resonance  and  the  resonance  in  the  "longi¬ 
tudinal"  system  are  superimposed  on  one  another  (Fig.  35).  The  amplitudes 
are  limited  by  nonlinear  damping  and  by  terms  of  higher 
order  (in  the  expression  for  the  nonlinear  function). 

The  factor  >/~3^  in  the  resonance  condition  stems 
from  the  nonlinearities  of  the  "longitudinal"  system.  In 
fact,  the  linear  coefficient  which  replaces  the  elasticity  of 
this  system  is  very  large  during  small  displacements  (in¬ 
finitely  large  during  displacements  tending  to  zero,  if  the 
rod  is  considered  incompressible)  and  approach  a  constant  value  c  when  the 
deflections  are  increasing.  This  process  is  shown  in  Fig.  36,  where  the 

reaction  of  the  system  to  the  longitudinal  displace¬ 
ment  is  plotted  along  the  vertical  axis. 

4.  We  shall  briefly  dwell  on  methods  for  preventing 
parametrically  excited  vibrations.  Besides  the  ob¬ 
vious  prevention  measure,  i.  e.  ,  the  reduction  or  the 
total  elimination  of  exciting  forces  and  also  the  re¬ 
moval  of  the  structure  from  the  dangerous  region 
by  meauis  of  changing  its  parameters,  the  following 
methods  are  recommended. 


Fig.  35 
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1)  Linear  damping.  This  prevention  measure  is  effective  only  when 
the  damping  is  large,  so  that  »  p.  For  >  p  parametric 
resonance  is  impossible.  For  a  choice  of  suitable  damping 
characteristics  one  should  use  the  equation  in  *9. 

2)  Nonlinear  damping.  The  introduction  of  a  "longitudinal"  damper 

reduces  the  amplitudes  of  the  parametric  vibrations  approxi¬ 
mately  inversely  proportional  to  This  method  should  be 

considered  only  as  a  helping  measure. 

3)  The  introduction  of  nonlinearities  by  springs  and  inertia  members. 
If  the  features  of  the  construction  allow  its  use,  such  a  method 
may  prove  to  be  useful.  One  must  remember,  however,  that 

the  influence  of  the  indicated  nonlinearities  is  opposite,  and  that 
there  exist  such  combinations  which  are  tinfavorable  to  the  de¬ 
crease  of  vibrations.  Consequently,  the  incre&se  of  a  nonlinear 
elasticity  of  the  system  does  not  always  reduce  the  amplitudes. 

The  same  is  true  for  nonlinear  inertia.  This  method  does  not 
permit  the  total  elimination  of  the  vibrations. 

•  22.  EXPERIMENTAL  VERIFICATION  OF  THE  THEORY 

1.  On  the  experimental  set-up  described  in  *6,  a  series  of  rods  were  in¬ 
vestigated  and  the  amplitudes  of  steady- state  vibrations  determined.  The 
data  of  the  experiment  were  compared  with  theoretical  results.  Some  of  these 
data  are  cited  below,  along  with  a  description  of  methods  for  the  determina¬ 
tion  of  the  nonlinear  characteristics. 

The  coefficient  of  nonlinear  elasticity  y  was  determined  from  the 
experiment  in  which  the  rod  was  loaded  with  a  longitudinal  force  exceeding 
the  critical  value.  In  *13,  the  formula 


/  = 


Ty 


(5.28) 


69 


was  derived.  Experiments  have  shown  that  Eq.  (5.  28) 
describes  quite  satisfactorily  the  dependence  of  the 
deflections  on  the  longitudinal  force  if  the  latter  is 
slightly  greater  than  the  Euler  buckling  force  (Fig.  37). 
From  Eq.  (5.28)  it  follows  that 


For  a  longitudinal  force  P  within  the  limits  F#  ^  F  ^  00065 

and  deflections  measured  in  the  middle  of  the  span,  we  obtained 
2  -4  2 

y/w  =  0.  42  X  10  cm  for  one  of  the  specimens. 

The  coefficient  of  nonlinear  damping  c  ^  was  determined  from  oscillo¬ 
grams  of  free  damped  vibrations. 

These  results  indicate  that  the  decrement  of  damping,  computed 
according  to  the  formula 


d  =  ^  In 

T  A{l  +  T) 


increases  with  amplitude.  The  experimental  values  for  the  above  mentioned 

specimen  are  plotted  in  Fig.  38.  The  increase  of  the  decrement  is  satis- 

_  2 

factorily  given  by  the  parabola  t  =  0.  10  +  0.  02  A  ,  from  which 

-2  - 1 

it  follows  that  €  =  0.  08  cm  sec 

The  coefficient  of  nonlinear  inertia  was  calculated  from 
the  formula 


3I*Mi 

4mP 


where  the  reduced  mass  is  found  from  elementary  considerations 

(Fig.  39). 
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We  will  introduce  the  following  notations:  let 
be  the  weight  of  the  specinnen,  be  the  weight  of 
loading  arm  CD  that  produces  the  constant  compo¬ 
nent  of  axial  force,  be  the  weight  of  the  moving 
support  of  the  vibrator  and  other  hardware,  and 
be  the  weight  of  the  additional  loading  on  the  end  of 
the  lever.  If  w  is  the  vertical  displacement  of  point 
B,  then  the  displacement  of  the  end  of  the  lever  arm 
will  be  w{a  +  b)/a.  The  axial  inertia  force  in  the  rod  is  then 


from  wliich 


Mv 


2^ 

t 


F,  4-lrCr, 


(“  +  »)*  Ip 

«•  r» 


For  the  case  where  Gj  =  6.  0  kg,  G^  =  31.0  kg,  Pj  =  30.  0  kg, 

P2  =  12.  0  kg,  a  =  40  cm,  b  =  130  cm,  f  =  170  cm,  we  obtain  Mj^g  =  4Z9.  0  kg. 
Consequently 


3.l4<-429 
*  4- 0.0  - 170* 


=  0.061  cm-*. 


The  coefficient  of  nonlinear  inertia  (more  accurately,  the  characteristic 
magnitude) 


j»  =  sn* 
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can  also  be  determined  from  the  experiment  on  the  damping  of  free  vibrations. 
According  to  the  approximate  formula,  Eq.  (4.  2Z),  the  frequency  of  the  non¬ 
linear  system  is 


Near  the  principal  resonance  n  «  1,  the  formula  takes  the  form 


0  =  0[l-{a*p) 


If  the  natural  frequencies  are  determined  for  various  amplitudes,  one  can 
compute  the  coefficient 


This  method,  however,  is  unreliable;  it  can  be  recommended  only  for 
systems  with  pronounced  nonlinearities. 

2.  The  resonance  curves  were  determined  in  such  a  manner  that  at  least 
five -sixths  of  the  experimental  points  are  within  the  limits  of  the  principal 
instability  region  (of  the  linear  approximation).  The  recording  was  made 
stepwise  for  increasing  and  decreasing  frequencies;  this  permitted  the  deter¬ 
mination  of  "overhanging"  vibrations. 

The  resonance  curve  for  one  of  the  specimens  is  given  in  Fig.  40. 

If  one  applies  the  equations  derived  from  theory,  one  must  consider  that  the 
amplitude  of  a  periodic  force  produced  by  the  vibrator  grows  proportionally 
to  the  square  of  the  exciting  frequency.  The  value  of  the  coefficient  at 
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n  =  1  is  denoted  by  ^  Then 

(5.29) 

This  expression  for  ^  should  be  substituted  in 
all  previously  derived  theoretical  formulas. 

We  will  give  the  results  of  the  calculations 
for  the  case  corresponding  to  Fig.  40  for  the 
following  values: 


Frequency  of  free  vibrations 

Euler  buckling  force 

Constant  component  of  the  axial  force 

Amplitude  of  the  variable  component  (at  n  =  1) 

—  ^t  54  2 

Excitation  parameter  p.  =  p-  y  =  jrggy  “  0.  212 

Reduced  longitudinal  mass 

Coefficient  of  nonlinear  inertia  X  =  ■  —  =  0.  044  cm  ^ 

(4)(6.0)(170)^ 


=  183  kg 
Pq  =  55  kg 
=  54.  2  kg 


Mj^g=  316  kg 


The  nonlinearities  of  the  elasticity  and  also  the  linear  and  nonlinear 
damping  are  so  small  that  they  can  be  neglected  in  the  calculation  of  ampli¬ 
tudes.  The  experimentally  determined  coefficient  of  linear  damping  obtained 
from  experiment,  for  example,  is  €  =  0.08  sec”^;  therefore,  the  critical 
coefficient  of  excitation  is 


73 
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The  coefficient  of  nonlinear  elasticity  y  is  small  compared  to  X12  >  sind 
the  influence  of  nonlinear  elasticity  can  therefore  also  be  neglected. 

The  boundaries  of  the  principal  region  of  excitation  are  determined 
from  Eq.  ( 1. 43), 


9. 


20 


Substitution  gives  61. 3  cps  and  75.  9  cps  for  the  lower  and  upper  boundaries, 
respectively.  Experimental  results  give  60.  8  cps  and  74.  2  cps  (Fig.  40). 
The  deviation  is  within  the  limits  of  the  accuracy  of  the  experimental  re¬ 
cording  of  the  oscillogram. 

The  vibration  amplitudes  are  determined  according  to  Eq.  (5.  22)  in 
conjunction  with  Eq.  (5.29): 


A*  = 


1  -„t(|  -^) 

Mil* 


2  — 

On  the  lower  boundary  of  the  region  of  excitation  n*  =  1/(1  +  p),  from 

2  —  ^ 
which  A  =  2p/X  follows.  In  our  case 


A 


2  •  0.212 
0.044 


0.84  cm . 


The  theoretical  results  are  plotted  in  Fig.  40  as  a  solid  line.  The 
somewhat  higher  values  for  the  experimental  amplitudes  near  the  upper 
boundaries  of  the  region  of  excitation  may  be  explained  by  the  fact  that  the 
vibrations  are  always  accompanied  by  more  or  less  intensive  beats.  How¬ 
ever,  in  the  figure  the  maximum  amplitudes  are  given  without  any  correction 
for  the  error  caused  by  the  beats. 
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Analogous  experiments  were  conducted  with  a  number  of  other  speci¬ 
mens,  where  the  constant  and  periodic  component  of  the  axial  force  were 
varied  for  each  of  them.  All  the  results  were  in  satisfactory  agreement 
with  the  theory.  The  experimental  data  are  shown  in  Figs.  41  and  42. 


( 


Fig.  42 


Figure  42  shows  the  variation  of  the  vibration  amplitudes  at  9  =  2(1  for  a 
changing  Pq.  In  the  case  of  increasing  Pq,  the  natural  frequency  (1  de¬ 
creases  and  the  variable  loading  P^  drops  accordingly. 
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